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ABSTRACT 


The  computer  program  presented  in  this  report  has  beep  designed 
specifically  for  the  analysis  of  the  blunt -trailing-edge  supersonic  com¬ 
pressor.  Beginning  with  flow  property  measurements  obtained  in  a  non¬ 
rotating  coordinate  system,  streamtube  boundaries  are  determined  at 
each  measuring  plane.  Then  mass-averaged  values  of  the  flow  proper¬ 
ties  in  each  streamtube  are  translated  to  a  coordinate  system  rotating 
with  the  compressor  rotor,  and  a  particular  one -dimensional  flow  model 
is  imposed  on  each  streamtube  to  describe  the  flow  process  through  the 
rotor.  The  flow  model  includes  analysis  of  shock  loss  and  sudden  expan¬ 
sion  loss  leading  to  an  estimate  of  the  additional  loss  occurring  within 
the  flow  field  of  the  rotor.  Various  other  calculations  are  presented 
which  are  aimed  at  the  analysis  of  data  for  accuracy  and  consistency. 
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j  Designation  of  the  axial  measuring  planes 
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=  2  Implies  the  measuring  plane  at  rotor  exit  in  which  traverse 

measurements  are  made  (This  counter  is  also  used  to  imply 
downstream  conditions  in  general. ) 

=  3  Implies  the  measuring  plane  at  rotor  exit  in  which  raJce 

measurements  are  made 

m  Streamtube  counter  implying  mass -averaged  conditions  of 

that  streamtube 

P  Profile 

PS  Passage  shock 

SE  Sudden  expansion 

SS  Suction  surface 

T  Based  on  stagnation  conditions 

t  Tip,  outer  casing 

6  Tangential  component 

SUPERSCRIPTS 

/ 

* 


Relative  to  the  rotor 
Critical  conditions 


Vlll 


AEDC-TR.69.42 


SECTION  I 
INTRODUCTION 


The  computer  program  described  herein  has  been  developed  for  the 
IBM  360/50  in  an  attempt  to  extend  the  capability  of  data  analysis  for 
the  supersonic  blunt-tr ailing- edge  compressor  at  the  Arnold  Engineer¬ 
ing  Development  Center.  An  example  of  a  typical  compressor  blade 
profile  is  shown  in  Fig.  1  (Appendix  I).  The  program  is  intended  pri¬ 
marily  to  give  some  insight  into  the  flow  process  through  a  compressor 
rotor  viewed  as  if  the  observer  were  rotating  with  the  angular  velocity 
of  the  rotor.  In  general,  theoretical  techniques  must  be  used  to  deter¬ 
mine  the  flow  in  a  rotating  passage  since  it  is  not  mechanically  possible 
to  place  a  probe,  within  the  rotor. 

The  problem  is  approached  basically  from  the  streamtube  point  of 
view.  If  flow  is  assumed  steady  and  axisymmetric  and  if  the  radial 
height  of  the  blade  is  divided  into  small  enough  increments  so  that  the 
radial  variation  within  the  increment  is  small,  then  it  shovild  be  possible 
to  consider  the  flow  through  the  rotor  as  one -dimensional  within  each 
streamtube.  Strictly  speaking,  this  method  demands  further  that  sec¬ 
ondary  flow  effects,  such  as  radial  mass  transfer,  be  negligible. 

The  model  for  the  composition  of  the  relative  total  pressure  losses 
in  a  streamtube  is  quite  similar  to  that  proposed  in  Ref.  1  for  flow  in 
transonic  compressors.  In  addition  to  the  losses  resulting  from  shock 
waves  and  the  viscous  effects  along  the  blade  profile  as  suggested  in 
Ref.  1,  the  blunt -trailing -edge  blading  may  be  expected  to  have  addi¬ 
tional  loss  caused  by  the  abrupt  area  increase  at  the  blade  trailing  edge. 
The  general  validity  of  such  a  model  has  been  demonstrated  in  Ref.  2. 

The  sections  of  this  report  include  a  description  of  the  computer 
input  (Appendix  n),  presentation  of  the  equations  used  in  the  computa¬ 
tion,  and  a  description  of  the  output  from  the  program  (Appendix  III) 

The  program  is  written  in  FORTRAN  IV. 


SECTION  II 

THE  EQUATIONS  AND  THE  COMPUTATIONAL  SCHEME 


In  general,  the  equations  of  this  section  outline  step-by-step  the 
process  of  the  computation.  The  procedure  is,  basically,  to  determine 
the  boundaries  and  centers  of  a  selected  number  of  steamtubes,  MM, 
using  data  measured  in  an  absolute  frame  of  reference,  then  to  mass 
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averaige  the  measured  data  between  streamtube  boundaries  and  to  consider 
these  average  values  as  concentrated  at  the  center  of  the  streamtube. 

The  averaged  values  are  then  translated  to  a  frame  of  reference  rotating 
with  the  angular  velocity  of  the  rotor.  The  flow  model  described  in  Refs. 

2  and  3  is  then  used  to  determine  additional  loss  unaccounted  for  in  the 
model.  Correlation  of  this  additional  loss  remains  the  basic  problem  of 
the  analyst. 

In  addition  to  the  main  objective  of  the  program,  described  above, 
several  additional  computations  have  been  included  to  provide  checks  on 
the  results  of  experiments  and  to  enlarge  the  use  of  the  existing  data. 


2.1  ASSUMPTIONS 

1.  The  flow  of  air  as  a  weightless,  thermally  perfect  gas 
is  assumed  to  be  frictionless,  time -steady,  adiabatic, 
and  axisymmetric  as  determined  from  discrete  point 
measurements  of  the  flow  variables. 

2.  Radial  mass  transfer  across  streamtube  boundaries  is 
ignored. 

3.  Five  discrete  points  are  selected  for  measurement  of 
the  flow  properties  in  front  of  the  wheel. 

4.  Behind  the  wheel  two  different  sets  of  measurements 
may  be  used  to  independently  describe  the  flow  at  the 
same  axial  plane  of  measurement.  In  the  tests  at  AEDC 
this  capability  has  been  used  to  evaluate  separate 
measurements  made  by  single-probe  traverse  and  m;jlti- 
probe  rake.  Traverse  measurements  may  be  supplied  at 
either  5  or  11  discrete  points  in  the  flow  field.  Rake 
measurements  are  supplied  at  only  5  points. 

5.  In  one  specific  situation  the  two  sets  of  measurements 
behind  the  wheel  may  be  obtained  at  different  axial 
planes.  This  particular  condition  is  differentiated  by 
the  capability  of  entering  exactly  13  discrete  point 
measurements  for  the  traverse.  The  set  of  data  for 
the  rake  will  still  contain  measurements  at  only  5 
discrete  points. 

6.  Discrete  point  measurements  of  total  pressure,  absolute 
flow  angle,  and  total  temperature  may  be  entered  in  all 
three  measuring  planes  at  5  points.  When  traverse  meas¬ 
urements  are  made  at  11  or  13  points,  only  total  pressure 
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and  absolute  flow  angle  are  obtained.  Total  temperature 
is  obtained  at  only  5  points  in  the  midportion  of  the  flow 
field  and  not  in  the  extra  points  traversed  in  the  boundary 
layer.  The  extra  total  temperatures  necessary  for  calcu¬ 
lations  are  obtained  in  the  program  by  extrapolation. 

7.  Static  pressure  is  measured  at  the  walls  at  the  axial  planes 
of  the  other  flow  property  measurements.  The  static 
pressure  is  assumed  to  vary  linearly  across  the  annulus  to 
obtain  values  at  the  discrete  points  of  measurement. 


2.2  DETERMINATION  OF  THE  PROFILES  FOR  STATIC  AND  TOTAL  PRESSURE, 

TOTAL  TEMPERATURE,  AND  THE  ABSOLUTE  FLOW  ANGLE 

The  inlet  measuring  plane,  designated  by  the  subscript  j  =  1,  is 
assumed  to  be  located  near  enough  to  the  leading  edge  of  the  rotor  blades 
so  that  flow  conditions  at  the  measuring  plane  may  be  considered  repre¬ 
sentative  of  conditions  immediately  ahead  of  the  rotor  leading  edge.  The 
downstream  measuring  planes,  designated  j  =  2  for  measurements 
obtained  by  traverse  and  j  =  3  for  measurements  obtained  by  rakes,  are 
assumed  to  be  far  enough  downstream  so  that  conditions  may  be  con¬ 
sidered  axisymmetric.  The  planes  of  measurement  are  assumed  to  be 
perpendicular  to  the  compressor  axis. 

In  each  measuring  plane  a  number  of  radial  positions,  Nj,  are  chosen 
as  points  at  which  measurements  are  made  of  total  pressure,  total  tem¬ 
perature,  and  absolute  flow  angle.  Static  pressure  is  measured  at  the 
wall  in  the  same  axial  plane  and  is  asstimed  to  vary  linearly  across  the 
annulus.  A  more  complete  description  of  the  measuring  techniques  and 
the  general  layout  of  the  supersonic  compressor  testing  facility  is 
described  in  Ref.  4.  For  j  =  1  or  3,  Nj  has  the  value  of  5.  For  3=2, 

Nj  may  have  the  value  of  5,  11,  or  13  corresponding  to  NN  in  the  input 
data.  The  radial  location  of  the  i-th  measuring  station  in  the  j-th  meas¬ 
uring  plane  is  given  by 

rji  =  rjh  +  Xji(rj,  -  rjh).  j  =  1,  2,  3,  i  =  1,2 . Nj.  0  <  xjj  <1  (1) 

where  rjh  is  the  radius  in  inches  of  the  inner  casing  wall  (nub  wall),  rjt 
is  the  radius  in  inches  of  the  outer  casing  wall  (tip  wall),  and  where  Xji 
may  be  interpreted  as  the  location  of  the  measuring  station  in  terms  of 
the  fractional  part  of  annulus  height.  The  quantities  rjh,  rjt,  and  xji 
aoce  input  data. 

Static  pressure  measurements  are  taken  on  the  inner  (h)  and  outer 
(t)  casing  walls.  The  assumption  of  straight  line  variation  across  the 
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axinvilus  implies  that  the  static  pressure  at  any  xji  may  be  calcijlated 
by  the  equation 

Pji  ~  Pjh  +  *ji  ^Pjl  ~  Pjhh  j  -  li  2,  3,  i  =  2 . ^  l,  0  <  Xji  <  1  ^2j 

When  Nj  equals  5,  the  static  pressure  is  calculated  for  two  additional 
points  at  Xji  equal  to  0.  05  and  0.  95. 


Because  of  the  physical  requirement  of  no  "slip"  at  a  wall,  the 
flow  velocity  at  a  wall  must  be  zero.  For  this  reason  the  static  and 
total  pressures  are  equivalent  at  a  wall;  i.  e. , 


Pji  -  Pjh, 

P  . 

=  Pit, 


i  =  h 
i  =  t 


j  =  1.  2,  3 


(3) 


The  drop  in  total  pressure  from  its  inviscid  value  occurs  rapidly  in  the 
wall  boundary  layer.  In  the  cases  of  N2  =  NN  =  11  or  13,  a  number  of 
total  pressure  measurements  have  been  made  in  the  boundary-layer 
region  so  that,  together  with  relations (3l[,  a  complete  total  pressure  pro¬ 
file  may  be  approximated.  In  the  cases  of  N-j  =  NN  =  5  the  data  are  in¬ 
sufficient  to  allow  accurate  mass  flow  calculations  across  the  annulus. 
Therefore,  two  additional  points,  located  at  xji  equal  to  0.05  and  0.  95, 
are  supplied  using  an  empirically  estimated  relation  based  on  the  more 
accurate  boundary -layer  surveys.  The  value  for  Nj  has  thus  been  in¬ 
creased  from  5  to  7  for  these  cases;  i.  e. ,  5  actual  measurements  in 
the  flow  field  entered  by  input  data  plus  2  extrapolated  points. 

In  every  case,  total  temperature  is  obtained  only  at  the  five  meas¬ 
uring  stations  in  the  midportion  of  the  annulus.  To  obtain  the  extra 
values  necessary  for  computation,  a  second  degree  polynomial  is  fit  to 
the  three  measured  points  nearest  the  locations  where  the  values  are 
needed. 

For  the  cases  when  it  is  necessary  to  estimate  the  total  pressure  at 
=  0.  05  or  0.  95,  it  is  also  necessary  to  supply  an  estimate  of  the 
solute  flow  angle  at  these  points.  The  boundary-layer  traverses  have 
sho^vn  that  it  is  sufficiently  accurate  to  assume  a  straight  line  fit  of  the 
two  points  nearest  the  annulus  wall  considered. 

2.3  DETERMINATION  OF  THE  ABSOLUTE  MACH  NUMBER,  TOTAL  AND  STATIC 
TEMPERATURE  CORRECTED  FOR  PROBE  RECOVERY  FACTOR,  AND  THE 
RATIO  OF  SPECIFIC  HEATS 

The  specific  heat  at  constant  pressure,  Cp,  is  slightly  temperature 
dependent.  For  air,  using  the  data  of  Ref.  5,  a  third-order  polynomial 
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was  found  to  fit  very  well  for  the  temperature  range  of  interest,  400  to 
900°R,  i.  e. , 

Cp  [Btu/lbm  °R1  =  a  -  bl  +  cd  -  dt“  (4) 


with 

a  =  0.241883 
b  =  -1.22830  X  10-* 
c  =  1.15098  X  10“‘ 

d  =  7.75163  X  10“"  (5) 


where  t  is  entered  in  degrees  Rankine.  Multiplying  the  result  by  g^J  to 
obtain  units  of  ft2/sec2°R  implies  that  the  ratio  of  specific  heats,  k, 
may  be  calculated  by 

K  =  Cp/(Cp  -  R)  (6) 

The  absolute  Mach  number  is  calculated  by 

Correction  of  the  measured  and  extrapolated  total  temperatures  for  the 
probe  recovery  factor,  RF,  may  be  accomplished  by 


T  (kM*  -  +  2) 

RF  (kM“  -  M*)  -  2 


(8) 


The  static  temperature  is  then  given  by 

.  -  T. 


(9) 


Beginning  with  an  initial  guess  on  the  static  temperature  in  Eq.  (4), 
an  iteration  process  is  csirried  out  with  Eqs,  (4)  to  (9)  until  Eq.  (9)  is 
satisfied  to  the  accuracy  of  the  machine  for  each  measuring  station,  i, 
in  every  measuring  plane,  j. 


2.4  DETERMINATION  OF  THE  TANGENTIAL  AND  AXIAL  COMPONENTS  OF  VELOCITY 
FOR  PLOTTING 


The  absolute  velocity  may  now  be  calculated  using  the  equation 


Vji  = 


Rt'ii  M 


ji  J 


=  1,  2,  3,  i  =  t,  1,  2, 


Ni.  h 


(10) 


The  axial  and  tangential  components  of  the  absolute  velocity  are  then 
given,  respectively,  by 


and 


Vaji  =  V ji  cos  oji 
V0ji  =  Vjisinoji, 


(11) 
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where  j  =  1,  2, 3,  i  =  t,  1,  2, . . .  ,Nj,  h. 

For  J  =  2  and  3  the  ratios  of  V^ji  and  Vgjj  to  the  maximum  axial 
velocity  are  plotted  versus  annulus  height  when  IPbOT  is  -Sr  or  3.  An 
example  of  the  plot  is  shown  in  Fig.  2.  HPT"  d 


2.5  DETERMINATION  OF  THE  RADIAL  POSITION  OF  STREAMTUBE  BOUNDARIES 
AND  CENTERS 


In  each  measuring  plane,  j,  the  annulus  is  to  be  divided  into  a 
number  of  areas,  MM,  through  which  equal  mass  flows,  where  MM  is 
input  data.  The  mass  average  of  the  flow  properties  in  each  of  the  MM 
areas  is  determined  and  is  considered  to  be  concentrated  at  the  radius 
representing  the  center  of  mass  of  each  area.  • 

The  axial  mass  flow  per  unit  area  is  given  by 

Pji  Mji  /,  ^ii~  1  11*2 

/T  - 2“  Mjijcosaji  (12) 

for  j  =  1,  2, 3,  i  =  t,  1,  2, . . .  ,Nj,  h.  The  total  mass  flow  at  each  meas¬ 
uring  plane,  mj,  is  cal  ciliated  using  the  equation 


The  quantities  mj  aire  printed  out  on  output  data  PAGE  5  in  Appendix  III 
for  flow  continuity  checks  and  comparison  with  other  methods  of  mass 
flow  measurements. 


Beginning  with  rji^it  =  rjt,  2(MM)  -  1  radii  are  determined  satisfying 
the  relation 


2(MM) 


(14) 


If  the  convention  is  adopted  that  i*  equals  zero  when  =  rjt,  then  the 
radii  of  the  streamtube  centers  are  given  when  i*  is  odd  and  the  stream - 
tube  boimdaries  are  given  when  i#  is  even.  The  radii  in  terms  of  annulus 
height  may  be  given  by 


- —,i  =  1,  2,  3,  i*  =  t,l,2 . 2(MM)  -  1,  h 

fjt  -  fjh 


(15) 


Both  rjisje  and  xji^  are  printed  on  the  output  data  PAGE  2  (Appendix  III). 


In  an  effort  to  check  the  accuracy  with  which  the  integrations  have  been 
performed,  the  integral  of  the  mass  flow  between  rj2(MM)-l  a-nd  rjh 
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is  calcvilated  and  compared  with  the  amount  which  shovild  remain;  i.  e. , 


2  (MM) 


I  Smj  = 


'•i2(MM)-  1 

2rr  J  r 

Nh 


(t) 


X  100,  j  =  1,2,3 


2(MM) 


(16) 


These  quantities  are  printed  on  output  data  PAGE  5  in  Appendix  III. 


2.6  DETERMINATION  OF  MASS-AVERAGED  STREAMTUBE  FLOW  PROPERTIES 

A  measured  flow  property  may  be  mass  averaged  between  the 
boundaries  of  a  streamtube  using  the  equation 


where  the  counter  iB  is  varied  so  that  radii  bounding  the  m-th  streamtube 
are  used.  The  quantity  q  may  be  the  static  pressure,  total  pressure,  or  T'o 
the  absolute  flow  angle;  and  q  represents  the  mass  average  of  that  quanti¬ 
ty.  The  results  of  mass  averaging  are  printed  on  output  data  PAGE  3 
in  Appendix  III.  The  equivalent  one -dimensional  average  over  the  annulus 
height  of  each  property  may  be  computed  by 

1  ““  _ 

"  "mm  £  ^  ~  1>  2,  3 

ni=  I 

These  results  are  also  printed  on  output  data  PAGE  3  in  Appendix  III. 

Using  the  mass-averaged  values  of  static  and  total  pressure  and 
total  temperature  and  Eqs.  (4),  (5),  (6),  and  (9),  an  iteration  cycle  is 
set  up  to  yield  an  estimate  of  the  mass-averaged  values  of  the  static 
temperature  and  the  ratio  of  specific  heats.  Then  Eq.  (7)  is  used  to 
calculate  the  streamtube  average  absolute  Mach  number. 

In  the  following  description,  a  quantity  with  a  subscript  m  implies 
the  streamtube  mass  average  of  that  quantity.  The  subscript  i  will  be 
retained  to  indicate  input  data  or  calculations  at  the  radial  measuring 
stations. 
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2.7  CALCULATIONS  TO  CHECK  THE  BALANCE  OF  THE  ENERGY  EQUATION  FOR 
ADIABATIC  FLOW 


If  the  total  enthalpy  is  defined 


H 


/  'pJT 


then  the  streamtube  total  enthalpy  may  be  calculated 


H 


jm 


=  aT 


jm 


(19) 


for  j  =  2,  3,  m  =  1,  2, ,  MM,  in  units  of  Btu/lbm,  where  the 

coefficients  a,  b,  c,  and  d  are  given  in  Eq.  (5).  The  change  in  total 
enthalpy  through  a  streamtube  is  the  given  by 


—  Hju,  —  Him,  j  —  2,  3,  m  =  1,  2,  ...  ,  MM  (20) 

These  quantities  are  printed  on  output  PAGE  5  in  Appendix  HI. 


The  streamtube  average  absolute  v^ocity  is  given  by  an  equation 
similar  to  (10)  with  tjj,  and  Mjj  replaced  by,  respectively,  the 
mass -averaged  values  of  the  ratio  of  specific  heats,  static  temperature, 
and  absolute  Mach  number.  The  streamtube  average  of  the  tangential 
velocity  is  then  given  ly 


for  j 


^  0jm 


V* 

*  jm 

irT~  “J" 


1,  2,  3,  m  =  1,  2, .... ,  MM,  in  units  of  Btu  sec/lbm-ft. 


(21) 


The  energy  equation  for  adiabatic  flow  with  work  addition  mav  be 
written 


d'H  =  nd'(rV0) 


where  the  symbol  d '  (  )  indicates  change  along  a  streamline, 
angular  speed,  fi,  in  units  of  sec”  is  given  by 

£2  =  - 

60  sec/min 


The  wheel 
(22) 


where  6  is  the  wheel  RPM  and  is  supplied  by  input  data.  The  degree  of 
balance  of  the  energy  equation  may  be  represented  by 

SEjm  =  Hjm  -  Him  -  7^^  (^j">  -  rim  Vgim^  ^23) 

^r  j  =  2,  3,  m  =  1,  2, . . . ,  MM.  The  percentage  of  balance  is  calculated 

(*  If), „■  fir  . ““  (24) 

This  representation  of  the  balance  of  the  energy  equation  is  printed  on 
output  data  PAGE  5  in  Appendix  III  along  with  the  change  in  total  enthalpy 
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through  the  corresponding  streamtube.  The  one -dimensional  average  of 
these  quantities  is  also  presented. 


2.8  CALCULATION  OF  STREAMTUBE  TOTAL  PRESSURE  RATIO  AND  EFFICIENCY 


The  streamtube  total  pressure  ratio  is  given  by  the  relation 
Pi 


R 


pjin  - 


jni 


1  rn 


,  j  =  2,  3,  m  =  I,  2 . MM 


(25) 


This  quantity,  along  with  its  equivalent  one-dimensional  average  calcu¬ 
lated  by  Eq.  (18),  is  printed  on  output  data  PAGE  6  in  Appendix  III. 


A  temperature -averaged  specific  heat  ratio,  « jniJ  is  calculated  by 


■  jm 


Cn  -  R 


^  dT 


IT.  T  1  m 


j  =  2,3,  m  =  1,2, 


MM 


(26) 


for  use  in  evaluating  streamtube  efficiency. 


Adiabatic  efficiency  is  defined  by  the  equation 


^  -sir 

where  6Hig  is  the  isentropic  enthalpy  rise  and  6H  is  the  actual  enthalpy 
rise.  Defining  the  isentropic  coefficient  Cjm  by 

K.  _  1 

i  =  2,  3,  m  =  1,  2,  ...  ,  MM  (27) 

then  the  isentropic  enthalpy  rise  is  given  by 


.  Ti„  [(Rpj„)N"  -  11  t  T|„  -  11 

♦  I  I'ln,  [(Bpin.)’'!"  -  11  -  i  T;_  KRpjJ*'..  -  11  (28) 

for  j  =  2,  3,  m  =  l,  2,...,  MM,  where  the  coefficients  a,  b,  c,  and  d 
are  given  in  Eq.  (5).  The  streamtube  adiabatic  efficiency  is  then  given 

by 

(5His)j„ 

»?adjm  =  - - - —  j  =  2,3,  m  =  1,2 . MM  (29) 

"jm  -  "Im 

This  quantity  and  its  equivalent  one -dimensional  average  are  printed  on 
output  data  PAGE  6  in  Appendix  III. 
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V 


2.9  DETERMINATION  OF  THE  FLOW  PROPERTIES  IN  RELATIVE  COORDINATES 


The  streamtube  average  tangential  component  of  the  relative  Mach 
number,  is  given  by 

2  ^  ^ 


M  0jlll  - 


M 


jm  Sin  fljni 


(30) 


for  j  =  l,  2,  3,  m  =  l,  2,...,  MM.  The  average  relative  Mach  number 
is  then  given  by 

M'jm  =  Wjm  cos  (31) 

for  the  same  range  of  j  and  m. 


The  streamtube  average  relative  flow  angle,  /3jjn.  is  given  by 

/  M 'a  \ 

iS 


jm  = 

arc  tan 

\Mjm  cos 

m  = 

1,  2. 

3,  .  .  .  , 

(32) 


The  relative  total  temperature  and  total  pressure  are  given,  respec¬ 
tively,  1:^ 


T  '  -  T  • 

^  jm  -  *  jm 


‘  ^  2  J 


^  jm 


3'  A  ''i'"  ^ 

jm  =  Pin'll  +  - 2 - 


(33) 

(34) 


for  j  =  1,  2,  3,  m  =  1,  2, . . . ,  MM. 


The  relative  Mach  number,  flow  angle,  total  pressure,  and  total 
temperature  are  printed  for  each  streamtube  and  each  measuring  plane 
on  output  data  PAGE  4  in  Appendix  III. 


2.10  CALCUUTION  OF  THE  OVERALL  RELATIVE  TOTAL  PRESSURE  LOSS 

In  the  preceding  sections  the  main  objective  has  been  to  locate  a 
number  of  streamtubes  and  determine  the  average  relative  flow  proper¬ 
ties  in  each  streamtube  before  and  after  the  rotor.  The  objective  be¬ 
comes  now  that  of  attempting  to  fit  the  flow  model  to  the  experimental 
results. 

To  the  benefit  of  centrifugal -type  turbomachinery,  a  flowing  fluid 
may  do  work  or  be  worked  upbn  by  simply  changing  the  distance  of 
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streamlines  from  the  axis  of  the  turbomachine.  To  correct  for  this 
fact  in  calculating  the  overall  relative  total  pressure  loss  for  an  axial- 
flow  machine,  the  ideal  (no  loss)  relative  total  pressure  ratio  is  given 


by 


(R'pjm) 


. . 


(35) 


6 

for  j  =  2,  3,  m  =  1,  2, . .  , ,  MM  from  Ref,  where  '<avjm 


*^jni  ~  m 


2 


where  MTjm  is  the  ratio  of  the  outlet  element  wheel  speed  to  the  inlet 
relative  stagnation  velocity;  i.  e. , 


M-, 


jn 


jrn 


j  =  2,  3  m  =  1,  2, 


MM 


I  cn 


(36) 


where  «Tm  i®  ibe  specific  heat  ratio  based  on  the  inlet  relative  total 
temperature.  The  overall  relative  total  pressure  loss  is  then  given  by 


(li 


jm 


=  (Rpjm)id 


(37) 


also  from  Ref.  6.  This  quantity  is  printed  on  output  data  PAGE  6  in 
Appendix  III. 


The  overall  relative  total  pressttre  loss  is  assumed  to  be  the  sum  of 
shock  loss,  profile  loss,  and  loss  attributable  to  sudden  area  expansion 
at  the  trailing  edge.  Shock  loss  and  sudden  area  expansion  loss  may  be 
estimated  under  certain  assumptions.  The  profile  loss  must  be  obtained 
by  subtraction  of  these  values  from  the  overall  loss,  and  then  correlated. 


2.11  CALCULATION  OF  SHOCK  LOSSES 

Pseudo-normal  shock  diffusion  is  inherent  in  the  conception  of  the 
blunt-trailing-edge  compressor  blading.  At  maximum  back  pressure 
the  losses  caused  by  a  pseudo-normal  shock  system,  neglecting  viscous 
interaction,  cannot  be  greater  than  the  loss  of  a  single  normal  shock  at 
the  same  inlet  conditions.  The  passage  inlet  Mach  number  of  the  com¬ 
pressor  is  assumed  to  be  the  average  of  the  blade  inlet  relative  Mach 
number,  Mijjj*  and  the  suction  surface  Mach  number  at  the  entrance  to 
the  passage,  MsSm-  The  Prandtl-Meyer  angle,  v,  of  Mssm  is  given 
by 

I^(MSS  m)  =  v(Mhn)  ^Im  -  ^Irr.  +  d>  SSm 
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if  M  >  1.  0,  where  i'(M  Prandtl-Meyer  angle  of  the  inlet 

relative  Mach  number,  where  amount  of  turning  re¬ 

quired  for  the  inlet  relative  flow  to  become  tangent  to  the  blade  surface 
at  the  leading  edge,  and  where  ^SSm  angle  of  expansion  from  the 

leading  edge  to  passage  entrance. 


The  Prandtl-Meyer  angle,  v,  for  a  Mach  number,  M,  is  calculated 

by 

V  =  ^  ■'  arctan  |  (M^  -  1)  -  arc  tan  -  l”  (39) 

Therefore,  is  easily  determined  by  this  equation;  but  ^'(Mgsjj^), 

given  by  Eq.  (38),  requires  iteration  to  obtain  Mgg^^^  from  Eq,  (39). 

The  inlet  blade  angle,  is  assumed  to  be  described  by  an 
arbitrary  function  of  the  form 

/S'l  =  A3  +  A^r  +  Ajt^  +  arctan  (A^t)  (40) 

r  r 

where  the  coefficients  are  input  data.  Suction  surface  expansion,  ^gg, 
is  assumed  to  be  described  by  a  similar  function 

<^cc  =  — ^  +  — —  +  Bj  +  B4r  +  Bjf^  (41) 

SO  f  p 

where  the  coefficients  are  input  data. 


The  passage  entrance  Mach  number  is  given  by 

M  Im  +  MSSm 

’’^avm  “  2 

This  quantity  is  printed  on  output  data  PAGE  6  in  Appendix  III. 


(42) 


The  total  pressure  ratio  across  a  normal  shock  with  upstream 
Mach  number  Mavm  is  given  by 

1  '<lm 

.  . .  1  '‘Im-l 

II  +  -I'"  -  1] 


R 


PSm 


■t 


2k 


'^Im  +  1 


'} 


'‘Im  -  j 

Mavm 

(Kin,  -  1)  M'avni  +  2, 

(43) 


The  loss  caused  by  the  pseudo-normal  shock  system,  neglecting 
viscous  interaction,  is  then  given  by 


PSm 


1  -  R 


P  Sm 


1  - 


Plr 


L  Plm  J 


,  m  =  1,  2,  ....  MM 


(44) 
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where  Cj  is  a  correction  factor  supplied  in  the  input  data  and  used  to 
adjust  for  the  fact  that  the  shock  system  may  not  produce  normal  shock 
loss.  Generally,  Ci  is  approximately  0.9. 

For  transonic  flow  with  high  subsonic  inlet  relative  Mach  numbers, 
shocks  may  also  be  expected  to  occur.  Total  pressure  losses  attrib¬ 
utable  to  these  shocks  may  be  estimated,  similar  to  Ref.  7,  by  calcu¬ 
lating  Mssm  based  on 

»'(MsSm)  =  0SSin 

rather  than  Eq.  (38).  If  Mavm.  calculated  by  Eq.  (42),  is  greater  than 
one,  shock  loss  is  calculated  through  Eqs.  (43)  and  (44);  if  is 

less  than  or  equal  to  one,  no  shock  loss  is  assumed  to  occur. 

The  passage  shock  loss,  psm#  printed  on  output  data  PAGE  6 
in  Appendix  III. 

The  static  pressure  through  the  normal  shock  system  is  given  by 

PPSm  =  c,p,„,  Tl  +  (M^vtn  -  1)]  m  =  1,2 _ MM  (45) 

L  '‘Im  +  1  J 

where  C2  is  a  correction  factor  supplied  in  the  input  data.  This  quan¬ 
tity  is  printed  on  output  data  PAGE  5  in  Appendix  III  for  comparison 
with  experiment  and  the  results  of  sudden  area  expansion  calculations. 


2.12  CALCULATION  OF  THE  LOSS  ATTRIBUTABLE  TO  SUDDEN  AREA  EXPANSION  FOR 

CYLINDRICAL  STREAM  SURFACES 

In  order  to  use  the  results  of  Ref.  8,  it  is  necessary  to  assume 
that  the  streamlines,  represented  by  the  streamtube  centers,  lie  on 
circular  cylinders  between  the  plane  of  the  trailing  edge  and  the  down¬ 
stream  measuring  planes  (j  =  2  or  3).  Such  an  assumption  implicitly 
demands  that  radial  shift  of  the  streamlines,  if  such  shift  occurs,  takes 
place  entirely  within  the  rotor. 

The  flow  is  assumed  to  leave  the  rotor  at  an  average  angle  equal  to 
the  rotor  exit  blade  angle,  |3g,  given  as  a  function  of  radius  by 

jSe  =  +  Dj  -r  D^r  +  Djr*  +  arctan  (D^r)  (46) 

r  r 

where  the  coefficients  are  input  data. 
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The  critical  relative  Mach  number  at  the  downstream  measuring 
planes  (j  =  2,  3)  is  given  by 


111  * 

•> 

1 

Jin 

1  - 

y 

"  jn 

(47) 


The  critical  Mach  number  of  the  relative  flow  at  the  trailing-edge  plane 
is  then  given  by 


M 


M:*’  ■  ,  i  =  2,  3,  ni  =  1,  2 . MM 


(48) 


where  ^  Ejm  calculated  for  r  =  the  radius  of  the  center  of  the 

m-th  streamtube  for  measurements  in  the  j-th  plane. 


The  relative  Mach  number  at  the  trailing-edge  plane  is  then  given 


by 


Mr- 

Kjni 


'<Ej-n  +  1 


,  j  =  2,  3,  ni  =  I,  2,  .  .  .  ,  MM 


(49) 


where  KEjm  calculated  by  Eqs.  (4),  (5),  and  (6)  for  some  initial  guess 
on  the  trailing-edge  static  temperature. 

/A2\ 

The  calculated  area  ratio  of  sudden  expansion,  >  where  A2  is 

the  streamtube  area  at  the  downstream  measuring  station  and  Ag  is  the 
flow  area  at  the  trailing  edge,  is  given  by 


'  ^  cjm 


Ejm  cos  ^  Ej,| 
jrn  COS  /3jni 


1  -h 

2  rj  J  m 


i  +  • 


Mj-m 


A  ,  K  •  vr* 

^  +  '^ajm  -’'jm 


■*"  ^aim  im  COS* 


-  COS  ^  f^jni  (50) 

for  j  =  2,  3,  m  =  1,  2,  ,  MM  where 

<^ajin  =  ('^Ejin  +  *^jin)/2 

The  relative  total  pressure  at  the  trailing  edge  is  given  by 

p^:j„  -  pj.  ::r  1  ,  -- 1  (sd 


M'j™ 

nos  jSjm 

1  ^ajm  1 

■M'l 

E  J  m 

'  +  2 

jm  Ejm 

COS  ^Ejni 

,  ajm  “  1 

1  1 

jm 

•  T 

_  2 
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The  static  pressure  at  the  trailing  edge,  based  on  the  sudden  area 
expansion  equations  for  cylindrical  stream  surfaces,  is  given  by 


P  Kjm 


M'jm  cos  /3jni 
M'Ejm  cos  ^Ejm 


1  M 


1 


'a  jm 


-  1 


MEio. 


(52) 


The  static  temperature  at  the  trailing  edge  is  then  given  by  an  equa¬ 
tion  similar  to  Eq.  (9)  involving  relative  flow  properties.  An  iteration 
cycle  is  set  up  for  the  calculations  of  Eqs.  (49)  through  (52)  until  the 
static  temperature  at  the  trailing  edge  converges  to  the  limit  of  the 
machine. 


The  calculated  area  ratio,  the  relative  total  pressure,  and  the 
static  pressure  at  the  trailing  edge  are  printed  on  output  data  PAGE  5 
in  Appendix  III, 


The  relative  total  pressure  loss  caused  by  sudden  area  expansion 
is  given  by 


'^SEjtn 


f*  Ejm  ~  P  jm 
P  1  m  —  p  Itn 


(53) 


These  quantities  are  printed  on  output  data  PAGE  6  in  Appendix  III. 


2.13  DETERMINATION  OF  THE  PROFILE  LOSS 


The  overall  relative  total  pressure  loss,  w',  of  Eq.  (37)  is  assumed 
to  be  the  sum  of  the  shock  loss,  sudden  area  expansion  loss,  and  profile 
loss  so  that  the  profile  loss  may  be  given  by 

<i>'pjin  =  s'jm  -  tu'PSm  -  W  SEjm  (54) 


This  may  be  put  into  a  more  standard  form  by  defining  the  profile  loss 
parameter,  PL, 


PL 


jm 


cup  cos  m 


2  a 


1  in 


(55) 


where  cr,  the  solidity,  is  assumed  to  have  a  distribution  as  a  function  of 
radius  given  by 

a  =  — L  +  — —  +  Fj  +  F«r  +  FsP  (56) 


where  the  coefficients 
is  calculated  based  on 
inlet  (j  =  1). 


are  input  data.  The  solidity,  for  use  in  Eq.  (55), 
the  radius  of  the  center  of  the  streamtube  at  the 
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The  results  of  the  Loss  calculations,  as  a  function  of  annulus 
height,  are  plotted  for  both  the  traverse  data  (j  =  2)  and  the  rake  data 
(j  =  3).  Examples  of  these  graphs  are  shown  in  Fig.  2.  The  points 
representing  shock  loss  are  plotted  at  the  annulus  height  of  the  stream- 
tube  centers  at  the  inlet  measuring  plane.  The  overall  loss  points  are 
plotted  at  the  annulus  height  of  the  streamtube  centers  for  outlet  meas¬ 
urements.  For  cylindrical  stream  surfaces,  the  points  locating  sudden 
expansion  loss  must  be  positioned  at  the  same  radial  location  as  the 
streamtube  centers  for  outlet  flow  measurements.  The  shaded  area 
represents  the  computed  profile  loss. 


2.14  CALCULATION  OF  THE  BLOCKAGE  FACTOR 


The  blockage  factor  is  defined  as  the  difference  between  the  geo¬ 
metrical  area  at  the  trailing  edge  through  which  flow  may  pass  and  the 
calculated  flow  area  divided  by  the  geometrical  area;  i.’  e. , 


Bj  =  - j  =  2.  3 


(57) 


The  geometrical  area,  Ag,  through  which  flow  may  pass  is  the  differ¬ 
ence  between  the  annulus  area  at  the  plane  of  the  trailing  edge  an'd  the 
area  occupied  by  the  blades.  Assuming  the  distribution  of  the  blade 
trailing-edge  thickness,  th,  may  be  written  in  the  form 

th  =  +  -j;-?  +  E,  +  E4r  Ejf’  (58) 

then  Ae  is  given  by 


Ae  =  (rr  -  -  r^Eh) 


rEt  -  i-Eh 
fEi fEh 


FEt  kf 

+  E2  In  7^  +  Ej(rEt  -  rth)  +  (r’Ei  - 

where  rg^.  and  r^j^  are,  respectively,  the  radii  of  the 
annulus  walls  at  the  plane  of  the  trailing  edge,  and  where  n  is  the  num¬ 
ber  of  blades.  The  coefficients  of  Eq.  (58),  the  radii  rEt  and  rg^.  and 
the  number  of  blades  are  all  input  data.  The  thickness,  th,  of  Eq.  (58) 
must  be  considered  the  arc  length  for  use  in  these  equations. 


r’Eh) 


(59) 


outer  and  inner 


The  streamtube  area,  A 2,  at  the  downstream  measuring  planes  is 
calculated  by. 


Ajjm  =  -  r=jiB+l) 


(60) 
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where  j  =  2,  3,  and  where  the  counter  iB  is  varied  so  that  the  stream- 
tube  boundaries  of  the  m-th  streamtube  are  selected. 


l^'rom  the  calculations  of  sudden  area  expansion  for  cylindrical 
stream  surfaces,  the  streamtube  area  in  the  plane  of  the  trailing  edge 
is  given  by 


(Ak)  = 

cj-n  (A,/Ae)^ 


(61) 


so  that  the  calculated  flow  area  at  the  trailing  edge  is  given  by 


MM 

(AK),j  =  Z  (Ak) 

m=  I 


cjm 


(62) 


for  j  =  2,  3.  The  blockage  factor  may  then  be  calculated  by  Eq.  (57), 
and  the  results  are  printed  on  output  data  PAGE  5  in  Appendix  III. 


2.1S  CALCULATION  OF  THE  DIFFUSION  FACTOR 


The  diffusion  factor,  D,  which  is  used  as  a  blade  loading  param¬ 
eter,  may  be  given  by 


—  1  — 


jtn 


M'l 


1  m 


/  ^ jm  ^jn 

W  ^Im  Ml 


2<^lni  M'lin 


(63) 


for  j  =  2,  3,  m  =  1,  2,  ,  MM.  The  results  of  this  calculation  are 

printed  on  output  data  PAGE  6  in  Appendix  III. 


SECTION  III 

DESCRIPTION  OF  INPUT  DATA 


The  input  data  necessary  for  the  operation  of  the  computer  program 
are,  basically,  an  analytical  description  of  the  rotor,  the  locations  of 
measuring  positions,  as  well  as  the  actual  measurements  taken.  The 
measurement  data  may  be  entered  in  one  of  two  ways:  (1)  listing  of  all 
measurements  or  (2)  through  the  use  of  data  arrays  produced  from  the 
data  reduction  program  for  the  supersonic  compressor  test  program. 

If  the  measurements  are  listed,  the  values  entered  must  represent  the 
time  average  values  of  the  measured  flow  properties  at  each  measuring 
station.  When  the  data  arrays  are  entered,  the  program  handles  the 
simple  averaging  process  where  necessary. 
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For  the  use  of  this  program  it  is  implicitly  assumed  that  outlet  rake 
measurements  and  traverse  measurements  involving  5  and  11  radial 
positions  are  located  at  the  same  plane.  Furthermore,  the  traverse  of 
13  radial  positions  may  be  located  in  an  axial  location  different  from 
the  plane  of  the  rake  measurements  (and  traverse  of  5  and  11  radial 
positions).  Sample  card  formats  for  input  data  are  given  in  Appendix  III. 


3.1  GROUP  DATA 


The  classification  of  group  data  includes  the  analytical  description 
of  the  rotor,  some  assumptions  on  the  behavior  of  the  flow  model,  and 
the  locations  of  the  measurement  stations  in  each  measuring  plane. 
This  data  must  be  entered  for  each  different  rotor  configuration,  but 
they  are  required  only  once  for  each  rotor  configuration. 


The  cards  are  numbered  below  with  the  computer  format;  the 
names  of  the  variables  are  given  along  with  their  description.  Every 
card  must  be  entered  in  the  order  given. 


IPRNT,  HEAD  1  FORMAT  (113,  19A4) 

=  0,  stop  computation;  therefore  a  blank  card  must  bb 
placed  at  the  end  of  all  input  data. 

=  1  for  complete  print  out  of  the  step  by  step  compu¬ 
tation. 

=  2  for  input  and  output  tables  only  (see  description  of 
.  Output  Data,  Section  IV). 

HEAD  1  An  identification  of  the  rotor  configuration.  There 

are  76  spaces  allowed. 


CARD  1. 
IPRNT 


CARD  2. 
MX 
IX 


RF 

NB 


MX,  IX,  RF,  NB  FORMAT  (213,  IE  12.4,  115) 

Number  of  sets  contained  in  the  group.  t 

=  0  if  a  static  pressure  tap  is  located  on  the  compres¬ 
sor  casing  (tip)  at  the  rotor  trailing  edge. 

=  1  if  no  static  pressure  tap  is  located  at  the  rotor 
trailing  edge. 

Recovery  factor  for  the  temperature  probes. 

Number  of  blades  on  the  rotor. 
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CARD  3. 
Kit 

Rih 

K2t 

R2h 

Rat 

Rsh 

CARD  4. 
■^i 


CARD  5. 
Bi 


Rit*  Rih.  R2t.  R2h*  Rat*  Rah 


FORMAT  (6E12.  4) 


Radius  of  the  compressor 
inlet  flow  measurements. 

Radius  of  the  compressor 
inlet  flow  measurements. 


casing  (tip)  at  the  plane  of 
(in. ) 

hub  wall  at  the  plane  of  the 
(in. ) 


Radius  of  the  compressor  casing  at  the  plane  of  the 
outlet  traverse  flow  measurements,  (in. ) 


Radius  of  the  compressor  hub  wall  at  the  plane  of  the 
outlet  traverse  flow  measurements,  (in.) 

Radius  of  the  compressor  casing  at  the  plane  of  the 
outlet  rake  flow  measurements.  This  value  must  be 
entered  only  if  R2t  i  Rst  R2h  ^  RShJ  i-  e. ,  only  if 
either  or  both  the  annulus  walls  converge  or  diverge, 
(in. ) 

Radius  of  the  compressor  hub  wall  at  the  plane  of  the 
outlet  rake  flow  measurements.  This  value  must  be 
entered  only  if  R2t  ^  Rst  R2h  ^  RShi  i-  e. ,  only  if 
either  or  both  the  annulus  walls  converge  or  diverge 
at  outlet,  (in.) 


Ai  FORMAT  (6E 12.  3) 

Coefficients  of  a  curve  describing  the  inlet  rotor  blade 
angle  P'l  in  radians  as  a  function  of  radius  from  the  axis 
in  the  form 

/S'  =  +  Aj  -  A^r  +  +  arctan  (Ajr) 

r  r 

If  any  of  the  coefficients  are  zero,  the  proper  12  spaces 
required  by  the  format  may  be  left  blank. 

Bj  FORMAT  (5E 12.  3) 

Coefficients  of  a  curve  describing  the  suction  surface 
expansion  ^33  as  a  function  of  radius  from  leading  edge 
to  passage  entrance  in  radians.  Coefficients  are 
entered  for  a  curve  of  the  form 

<SsS  =  -^  +  -^  +  B3  +  B,r  .  B,r* 
r  r 

If  any  of  the  coefficients  are  zero,  the  proper  12  spaces 
required  by  the  format  may  be  left  blank. 
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CARD  6.  Cl,  C2  FORMAT  {2E 12.  3) 

Cl  Passage  shock  loss  adjustment  coefficient.  (0.9  is 

suggested.  See  Eq.  (44)) 

C2  Normal  shock  static  pressure  rise  adjustment  coeffi¬ 

cient.  (See  Eq.  (45)) 


CARD  7. 
Di 


Di  FORMAT  (6E  12.  3) 

Coefficients  to  describe  the  exit  blade  angle  in 
radians  as  a  function  of  radius.  The  assumed  form  is 

^  E  =  +  Dj  -  D^r  +  Dj  F  +  arctan  (Dji) 

If  any  of  the  coefficients  are  zero,  the  proper 
12  spaces  required  by  the  format  may  be  left  blank. 


CARD  8.  Ei  FORMAT  (5E12.  3) 

El  Coefficients  of  the  curve  to  describe  the  blade  trailing 

edge  thickness,  th,  as  a  function  of  radius.  Coefficients 
are  entered  for  a  curve  of  the  form 

K, 

ih  =  — —  -i-  —  -H  E,  E^r  4^  EjT^ 
r  r 

If  any  of  the  coefficients  are  zero,  the  proper 
12  spaces  required  by  the  format  may  be  left  blank. 


CARD  9. 
Fi 


Fi  FORMAT  (5E 12.  3) 

Coefficients  of  the  curve  to  describe  the  blade 
solidity,  o,  as  a  function  of  radius  in  the  form 

<7  =  —  +  —  -Fj  +  F*r  +  Fj  r' 
r*  r 

If  any  of  the  coefficients  are  zero,  the  proper 
12  spaces  required  by  the  format  may  be  left  blank. 


CARD  10.  Xii  FORMAT  (5E 12.  3) 

Xii  Fractional  part  of  annulus  height  locating  positions  of 

inlet  flow  measurements.  Listing  is  made  from 
casing  to  hub,  i.  e.,  decreasing  annulus  height.  Only 
5  are  entered. 
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CARD  11. 
X52i 

CARD  12, 

Xll2i 

CARD  15, 
X132i 

CARD  18. 
Xsi 

CARD  19. 

Zi 


CARD  20. 
R-ti 


X52i  FORMAT  (5E12.  3) 

Fractional  part  of  annulus  height  locating  traverse 
flow  measurements  for  outlet  when  5  positions  are 
recorded.  Listing  is  entered  from  casing  to  hub,  i.  e., 
decreasing  annulus  height. 

13,  14.  Xll2i  FORMAT  (5E12.  3) 

Fractional  part  of  annulus  height  locating  traverse 
flow  measurements  for  outlet  when  11  positions  are 
recorded.  Listing  is  entered  from  casing  to  hub,  i.  e. , 
decreasing  annulus  height. 

16,  17.  X132i  FORMAT  (5E12.  3) 

Fractional  part  of  annulus  height  locating  traverse 
flow  measurements  for  outlet  when  13  positions  are 
recorded.  Listing  is  entered  from  casing  to  hub, 
i.  e. ,  decreasing  annulus  height. 

Xgi  FORMAT  (5E12.  3) 

Fractional  part  of  annulus  height  locating  rake  flow 
measurements  for  outlet.  Only  5  are  entered.  List¬ 
ing  is  made  from  casing  to  hub,  i.  e. ,  decreasing 
annulus  height. 

Zi  FORMAT  (7E 11,  4) 

Axial  position  of  the  trailing  edge  of  the  rotor  and  the 
wall  static  pressure  taps  between  the  trailing  edge 
and  the  plane  of  the  downstream  rake  measurements, 
inclusive.  The  listing  is  entered  in  successive  order. 
A  maximum  of  7  may  be  entered.  If  there  are  less 
than  seven  values,  spaces  are  filled  until  the  required 
number  has  been  entered  and  blank  spaces  are  left  at 
the  end  of  the  card.  (in. ) 

Rti  FORMAT  (7E1 1.4) 

Radii  of  the  casing  walls  at  each  axial  position  cor¬ 
responding  to  the  static  taps  as  mentioned  on 
CARD  19.  This  card  is  entered  blank  if  Rti  is  con¬ 
stant  since  this  value  has  been  entered  on  CARD  3. 

(in. ) 
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CARD  21.  Rhi  FORMAT  (7E 11. 4) 

Rjji  Radii  of  the  hub  wall  at  each  axial  position  correspond¬ 

ing  to  the  static  taps  as  mentioned  on  CARD  19.  This 
card  is  entered  blank  if  Rh  is  constant  since  this  value 
has  been  entered  on  CARD  3.  (in. ) 

3.2  SET  DATA 

One  card  is  required  to  give  the  number  of  subsets  included  in  the 
set  and  to  allow  naming  of  the  set.  A  set  is  usually  interpreted  as  the 
data  required  to  produce  a  constant  rotor-speed  operation  character¬ 
istic.  The  card  is  placed  before  the  subset  data  included  in  that  set. 

CARD  22.  ISTART,  HEAD  2  FORMAT  (113,  19A4) 

ISTART  Number  of  subsets  included  in  the  set. 

HEAD  2  An  identification  of  the  set.  There  are  76  spaces 

allowed. 

3.3  SUBSET  DATA 

This  classification  of  data  includes  the  actual  flow  measurements. 
The  data  may  be  entered  in  alternative  methods  depending  on  whether 
data  are  to  be  listed  in  full  or  advantage  is  taken  of  the  data  arrays  out¬ 
put  by  a  data  reduction  program  for  the  supersonic  compressor  test 
facility.  A  maximum  of  12  subsets  may  be  included  with  each  set. 

CARD  23.  HEAD  3  FORMAT  (20A4) 

HEAD  3  An  identification  of  the  subset.  There  are  80  spaces 

allowed. 

CARD  24.  NN,  MM,  IPT,  INP,  RPM  FORMAT  (413,  IE  12.4) 

NN  Number  of  measuring  stations  at  the  outlet  traverse 

measuring  plane. 

MM  Number  of  streamtubes  sought;  there  may  be  a  maxi¬ 

mum  of  10. 

IPT  =  0  if  no  plots  are  desired. 

=  1  if  only  pressure  loss  plots  are  desired. 

=  2  if  only  velocity  ratio  plots  are  desired. 

=  3  if  both  velocity  ratio  and  pressure  loss  plots  are 
desired. 
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CARD  24.  (Continued) 

INP  =  0  if  input  is  by  data  array  translation;  CARDS  41 

through  46  are  expected  to  be  supplied  for  each  subset. 

=  1  if  input  is  by  a  full  listing  of  averaged  data; 

CARDS  25  through  40  are  expected  for  each  subset. 

RPM  Wheel  rotations  per  minute;  entered  only  if  INP  =  1. 


CARDS  25  through  40  are  expected  if  INP  =  1. 


CARD  25. 
Pit 

Plh 

P2t 

P2h 

Pst 

Psh 


Pit'  Plh'  P2t'  P2h’  PSt'  PSh  FORMAT  (6E 12.  4) 

Static  pressure  at  outer  casing  at  the  inlet  measuring 
plane.  (Ibf/in.^) 

Static  pressure  at  inner  casing  at  the  inlet  measuring 
plane.  (Ibf/in.^) 

Static  pressure  at  outer  casing  for  the  outlet  traverse 
measuring  plane.  (Ibf/in.^) 

Static  pressure  at  inner  casing  for  the  outlet  traverse 
measuring  plane.  (Ibf/in.^) 

Static  pressure  at  outer  casing  for  the  outlet  rake 
measuring  plane.  If  the  outlet  traverse  and  rake 
measurements  are  made  at  the  same  plane,  this  value 
is  not  entered.  (Iflf/  in.  2) 

Static  pressure  at  inner  casing  for  the  outlet  rake 
measuring  plane.  If  the  outlet  traverse  and  rake 
measurements  are  made  at  the  same  plane,  this 
value  is  not  entered,  (lbf/in.2) 


CARD  26.  p^i  FORMAT  (7E11. 5> 

Outer  casing  static  pressure  entered  in  a  manner  cor¬ 
responding  to  the  Z's  of  GROUP  DATA  CARD  19. 
Leave  blanks  for  unmeasured  values.  (Ibf/in.^) 


CARD  27. 

Phi 

CARD  28. 
Pli 


Phi 


FORMAT  (7E11.  5) 

Inner  casing  static  pressures  entered  in  a  manner  cor¬ 
responding  to  the  Z*s  of  GROUP  DATA  CARD  19. 
Leave  blanks  for  unmeasured  values.  (Ibf/in.^) 


li 


FORMAT  (5E12.  3) 


Inlet  total  pressure.  (Ibf/in.^) 
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CARD  29  (30,  31)  P2i  FORMAT  (5E 12.  3) 

P2i  Outlet  traverse  total  pressure.  If  NN  =  5,  one  card 

is  sufficient;  if  NN  =  11  or  13,  three  cards  are  re¬ 
quired.  (lbf/in.2) 


CARD  32. 

Psi 


P 


3i 


Outlet  rake  total  pressure. 


FORMAT  (5E12.  3) 

(lbf/in.2) 


CARD  33. 

Tli 


Tii  FORMAT  (5E 12.  3) 

Inlet  total  temperature.  (®F) 


CARD  34. 
T2i 


T2i  FORMAT  (5E12.  3) 

Outlet  traverse  total  temperature.  (®F) 


CARD  35. 

Tsi 


CARD  36. 

«li 


T3i  FORMAT  (5E12.  3) 

Outlet  rake  total  temperature.  (“F) 

aji  FORMAT  {5E 12.  3) 

Inlet  absolute  flow  angle  in  degrees. 


CARD  37  (38,  39)  Qr2i  FORMAT  (5E12.  3) 

0'2i  Outlet  traverse  absolute  flow  angle  in  degrees.  If 

NN  =  5,  one  card  is  sufficient;  if  NN  =  11  or  13, 
three  cards  are  necessary. 


CARD  40.  a^i  FORMAT  (5E 12. 3) 

ffsi  Outlet  rake  absolute  flow  angle  in  degrees.  This  card 

must  be  supplied  only  if  NN  =  13  (CARD  24), 

CARDS  41  through  46  are  expected  if  INP  of  CARD  24  equals  zero. 


CARD  41. 
‘^Si 


CARD  42, 

Pti 


a^i  FORMAT  {5E12.  3) 

Absolute  flow  angles  in  degrees  for  the  rake  measure¬ 
ments.  This  card  must  be  entered  only  if  NN  =  13; 
i.  e. ,  if  NN  =  13  (CARD  24). 


Pti  FORMAT  {5E 12.  3) 

Traverse  total  pressures  near  the  outer  casing.  This 
card  is  necessary  only  if  NN  =  11  or  13.  Begin  listing 
with  the  measurement  nearest  the  outer  casing. 
(lbf/in.2) 
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CARD  43.  Phi  FORMAT  (5E 12.  3) 

Phi  Traverse  total  pressures  near  inner  casing.  This 

card  is  necessary  only  if  NN  =  11  or  13.  Begin  list¬ 
ing  with  the  measurement  nearest  the  mainstream 
values,  (lbf/in.2) 


CARD  44.  ati 


FORMAT  (5E12.  3) 


ttti  Traverse  absolute  flow  angles  in  degrees  nearest 

outer  casing.  This  card  is  necessary  only  if  NN  =11 
or  13.  Begin  listing  with  measurements  nearest  the 
outer  casing. 


CARD  45.  Qfhi 


FORMAT  (5E12.  3) 


tthi  Traverse  absolute  flow  angle  in  degrees  near  inner 

casing.  This  card  is  necessary  only  if  NN  =  11  or  13. 
Begin  listing  with  measurement  nearest  the  main¬ 
stream  values. 


Arrangement  of  the  measured  values  for  CARDS  42  through  45  are 
entered  in  such  a  manner  that  the  mainstream  measurements  in  the 
data  arrays  may  be  placed  in  between  these  values. 

CARD  46.  NX,  COR(NX),  M  FORMAT  (5(  113,  1E12.  3),  115) 

NX  Channel  number  to  be  corrected. 

COR(NX)  Correction  to  be  substituted  in  channel  NX. 

M  =  0  if  5  corrections  are  sufficient;  i.  e. ,  one  card. 

=  1  if  another  card  of  corrections  is  necessary. 

CARD  46  may  be  entered  whenever  data  arrays  are  supplied.  This 
allows  substitution  for  incorrect  values  in  the  data  arrays.  The  card 
must  be  entered  blank  if  there  are  no  corrections. 

The  complete  program  is  actually  divided  into  two  main  sections; 
the  first  simply  prepares  the  data  for  entry  into  the  main  calculations 
described  in  Section  II.  In  the  process  the  prepared  data  are  stored  on 
tape  and  there  exists  the  option  of  recalling  this  tape  to  run  the  second 
part  of  the  program  independent  of  the  first  to  obtain  a  different  output. 
This  option  is  controlled  by  the  following  two  cards. 


25 


AEDC.TR.69.42 


CARD  47. 

ITAPE 

FORMAT  (113) 

ITAPE 

=  1  if  input  to  second  part  is  by  cards. 
=  2  if  input  to  second  part  is  by  tape. 

CARD  48. 

IPRNT 

FORMAT ,( 113) 

IPRNT 

=  0  to  stop  computation. 

=  1  for  complete  printout  of  step  by  step  computation. 

=  2  for  input  and  output  tables  only. 

This  option  is  useful  in  tracing  errors.  CARD  48  supercedes  the 
IPRNT  of  CARD  1.  These  two  cards  are  entered  once,  and  a  blank  card 
follows  to  indicate  program  completion. 

SECTION  IV 

DESCRIPTION  OF  THE  OUTPUT  DATA 


The  output  of  this  program  may  include 

1.  a.  Results  of  step  by  step  calculations,  or 
b.  Six  pages  summarizing  the  main  results. 

2.  Graphs  prepared  on  a  Calcomp  plotter  for: 

a.  The  ratio  of  the  axial  and  tangential  components  of  the  ab¬ 
solute  velocity  downstream  of  the  rotor  to  the  maximum 
axial  velocity  in  function  of  annulus  height,  or 

b.  The  composition  of  the  overall  relative  total  pressure  loss 
according  to  the  flow  model  in  function  of  annulus  height, 
or  both. 

In  the  following  pages  a  brief  description  of  Item  1.  b.  will  be  pre¬ 
sented.  An  example  of  this  output  is  shown  in  Appendix  III.  The  results 
of  the  step  by  step  calculations  were  used  for  checkout  purposes  and 
continue  to  be  helpful  in  locating  errors  in  input  data;  however,  detailed 
description  is  unnecessary.  Examples  of  the  graphs  are  shown  in  Fig.  2. 

Output  data  PAGE  1  presents  a  brief  summary  of  the  program  input 
data. 

Output  data  PAGE  2  presents  the  calculated  boundaries  and  centers 
of  the  MM  streamtubes  for  each  measuring  plane.  These  are  presented 
both  as  the  fractional  part  of  annulus  height  and  the  radius  in  inches 
from  the  compressor  axis. 
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Output  data  PAGE  3  presents  a  summary  of  the  mass -averaged 
values  of  the  streamtube  flow  variables  in  the  absolute  coordinate  sys¬ 
tem  for  each  of  the  measuring  planes.  Streamtubes  are  numbered  from 
tip  to  hub  in  the  first  column.  The  second  column  repeats  the  stream- 
tube  center  in  the  fractional  part  of  the  annulus  height.  The  mass- 
averaged  values  of  total  pressure  in  units  of  lbf/in.2  are  presented  in 
the  third  column.  Similarly,  the  fourth  column  presents  static  pressure 
in  units  of  Ibf/in.  the  fifth  column,  total  temperature,  in  “R;  the  sixth 
column,  absolute  flow  angle,  in  degrees.  The  one -dimensional  average 
in  each  measuring  plane  is  also  presented  for  the  flow  variables. 

Output  data  PAGE  4  presents  some  of  the  main  calculated  flow  con¬ 
ditions  and  their  one -dimensional  averages.  The  first  and  second 
columns  again  present  streamtube  number  and  the  location  of  the 
streamtube  center  in  the  fractional  part  of  annulus  height.  The  third 
and  fourth  columns  present  the  streamtube  absolute  Mach  number  and 
the  relative  Mach  number,  respectively.  The  fifth  column  presents 
the  streamtube  absolute  total  enthalpy  in  units  of  Btu/lbm.  The  sixth 
column  presents  the  relative  flow  angle  at  each  of  the  measuring  planes 
in  degrees.  The  seventh  column  presents  the  blade  angle  in  degrees 
calculated  at  the  radius  of  the  streamtube  center.  The  eighth  column 
presents  the  relative  total  pressure  in  lbf/in.2;  the  ninth  column,  the 
relative  total  temperature  in  “R. 

Output  data  PAGE  5  presents  the  total  integrated  mass  flow  in 
Ibm/sec  at  each  measuring  plane  along  with  the  estimated  percentage 
error  in  determining  streamtube  boundary  and  center  location.  It  pre¬ 
sents  the  calculations  exhibiting  the  balance  of  the  energy  equation, 
various  calculated  values  of  static  and  total  pressure  at  the  rotor  trail¬ 
ing  edge,  the  results  of  the  blockage  factor  calculation,  and  the  stream- 
tube  area  ratio  from  the  sudden  expansion  calculations. 

Output  data  PAGE  6  presents  the  principal  results  of  the  stream- 
tube  calculations  and  generally  involves  the  use  of  measured  data  up¬ 
stream  and  downstream  of  the  wheel.  In  addition  to  the  streamtube 
number  as  a  reference,  presented  are  the  streamtube  average  and  the 
one -dimensional  average  of  the  absolute  total  pressure  ratio,  adiabatic 
efficiency,  the  overall  relative  total  pressure  loss  coefficient,  the 
passage  entrance  Mach  number,  the  passage  shock  loss  based  on  this 
Mach  number,  the  loss  in  relative  total  pressure  attributable  to  sudden 
expansion,  the  calculated  (arithmetically)  profile  loss,  the  profile  loss 
parameter,  and  the  diffusion  factor. 

Computation  time  for  one  case  is  approximately  two  minutes. 
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APPENDIX  III 
OUTPUT  DATA 
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0.0  9.8435  5.0  9.9000  0.0  9.9000 


V 
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FROM  RUN  CA  11  0004  SFLECTFO  €X*MPLES  FOR  TEST  CASES 

DATA  GROUP  1  l.ON  MIN  EXHIBITING  OUTPUT  FOR  5  STREAMTUBES  WHEN  NN«5 

SUMMARY  OF  MASS  AVERAGEO  STREAMTUBE  FLOM  VARIABLES 

PAGE  3 

RADIAL 

RATIO  TO 

TOTAL 

STATIC 

TOTAL 

ABSOlUTF 

POSITION 

PASSAGE 

PRESSURE 

PRESSURE 

TEMP 

FLOW  ANGLE 

HEIGHT 

INLFT 

1 

r*.  ARPrt 

14.7865 

11.7029 

528.2647 

0.0 

2 

''.6900 

14.0434 

11.T418 

526.7777 

n,n 

3 

D.B14B 

14.9375 

11.TB19 

526.2191 

0.0 

4 

D.123T 

14.9228 

11.8231 

526.9851 

0.0 

5 

0.12T1 

14.8740 

11.8653 

529.D062 

0.0 

AVERAGE 

14. 8928 

11.7830 

527.4506 

p.p 

OUTLET 

TRAVERSE 

1 

O.7O6O 

18.8101 

12.7366 

680.7538 

45.5406 

2 

0.«.1P?I 

21.8457 

12.6960 

696.1969 

34.6380 

3 

0.4S7O 

23.046B 

12.6561 

701.4062 

34.6209 

- 

4 

P.1032 

21.2484 

12.6146 

689.9601 

15.3193 

8 

0.1223 

19.T691 

12.5701 

678.4390 

34.2794 

AVERAGE 

20.9440 

12.6547 

689.3512 

16.8798 

rDUTLPT 

RAKE 

1 

O.TSS'i 

19.2960 

12.7330 

684.2B09 

43.FS80 

2 

O.BIST 

22.1139 

12.6963 

691.9601 

34.8849 

3 

0.4frT6 

22.7793 

12.6561 

693.1563 

34.6406 

..  .  . 

4 

0.3D21 

21.0746 

12.6144 

689.0272 

35.3176 

5 

0.1214 

19.7785 

12.5698 

679.1307 

34.2702 

AVE<<AGE  2l#P08A  36*5?84 


AEDC-TR.69.42 


FROM  RUN  C»  31  0094  SFLECTEO  EXAMPLES  FOR  TEST  CASES 
DATA  GROUP  1  l.ON  MIN  EXHIBITING  OUTPUT  FOR  5  STREAMTUBES  NHEN  NN»5 

calculated  flow- conditions  ” 


03 


RADIAL 

RATIO  TO 

ABSOLUTE 

RELATIVE 

TOTAL 

RELATIVE 

BLADE 

RELATIVE 

RELATIVE 

POSITION 

PASSAGE 

MACH 

NACH 

ENTHALPY 

FION 

ANGLE 

TOTAL 

TOTAL 

HEIGHT 

NUMBER 

NUMBER 

ANGLE 

PRESSURE 

TEMPERATURE 

INLET 

1 

0.9820 

0.3877 

1.5807 

126.7801 

68.1746 

61.0340 

48.3568 

741.4553 

2 

0.6999 

0.5970 

1.5607 

126.4237 

67.5112 

60.5495 

47.1197 

731.6788 

..  .3 _ 

.QjS146 

0.5921 

1.5317  . 

024.2897 

67.2583 

60.0406 

722.8362 

4 

0.3237 

0.5863 

1.4997 

126.4734 

66.9872 

59.4998 

43.3998 

715.2469 

5  - 

(U127I  ^ 

0.5774 

1.4640 

126.9579 

66.7729 

58.9241 

41.3657 

70«.fl473 

AVERAGE 

0.5881 

1.5274 

126.5850 

67.3408 

60.0096 

45.1079 

724.0129 

outlet 


TRAVERSE 

.1 

0.7969 

0 17679 

0,9341 

163.4420 

54.8433 

30.7815 

22.3548 

715.0986 

2 

0.6193 

0.9160 

1.0850 

167.1701 

46.0029 

30.3154 

26.6014 

3 

0.4679 

0.9666 

1.0842 

168.4284 

42.8055 

29.9147 

26.4946 

729,8666 

A 

0.3032 

0.8965 

1.0414 

165.6640 

45.3831 

29.4751 

25.0703 

723.2945 

5 

0.1223 

0.8312 

1.0280 

162.8834 

48.0754 

28.9877 

24,5838 

721.9590 

AVERAGE 

0.B757 

1,0345 

165.5176 

47.4220 

2P.B9AP 

25.0210 

725.3296 

OUTLET 

RAKE 

i_  _ 

0.7855 

0.7944 

0.9568 

164.2931 

53.0107 

30.7516 

22.9256 

718.7466 

2 

0.6167 

0.9270 

1.0876 

166.1469 

45.6151 

30.3086 

26.6877 

7’?0.0fl80 

3 

0.4676 

0.9563 

1.0854 

166.4357 

43.5426 

29.9139 

26.5347 

724.0073 

A 

0.3023 

0.5889 

1.0397 

165.4388 

45.7681 

29.4726 

25.01  86 

723.5840 

5 

0.1214 

0.8317 

1.0276 

163.09B6 

48.018A 

28.9854 

24.5709 

7?2.7000 

AVERAGE 

0.8797 

1.0394 

165.0826 

47.1910 

29.8864 

25.1475 

723.8252 

PAGE  4 


CP-69-ai-Da3V 


COMMENTS 


o»GE 


MASS  FLOW _ 

PERCENT  Mass  LOSS  IN  DETERMINING  RADIAL  POSITIONS 

PERCENT  ENERGY  UNBALANCE 


STREAMTURF 


OUTLET  TRAV 


INLET  OUTLET  TRAV. 

gP.RTfeP  DO _ ?1.6«An  00 

<>.6530-09  2.0670-07 


OUTLET  RAKE 


OUTLET  RAKE 
21.5620  op 


OELTAIH) 

DELTAIEI/OELTAIHI 

OFLTAJHI 

DELTAIEI/OELTAIHI 

1 

36.662 

-1.452520 

01 

37.513 

-1.15'»22n 

"1 

2 

40.746 

4.916170 

00 

39.723 

1.25*510 

op 

3 

42.139 

4.984550 

00 

40.146 

1. th^orO 

fyn 

4 

39.191 

5.018740 

00 

?R*965 

S.?47rR0 

fsiy 

5 

35.926 

8.127710 

36.141 

R,«07Q5n 

r>n 

AVERAGE  38.933 

1.704400 

?0 

^e.49R 

1.054260 

r\fy 

STATIC  PRESSURE 

AT  TRAILING  EDGE 

RELATIVE 

TOTAl  PRFSSURE 

AT 

TRAILING  1 

FROM 

NORMAL  SHOCK  EQUATIONS  FROM  SUDDEN  EXPANSION 

EQUATIONS  FROM  SUDDEN  EXPANSION 

EQUATIONS 

streamtube 

PAISAGF  ENTRANCE 

OUTLET  TRAV. 

OUTLET  RAKE 

OUTLET  TRAV. 

OUTLFT  RAKF 

1 

45. IDO 

5.596 

5.456 

29.710 

2 

44.017 

4.904 

4.997 

27.237 

27.718 

3 

42.B77 

5.815 

5.501 

26.770 

26.872 

4 

41.623 

5.115 

5.014 

26.240 

26.303 

5 

40.280 

4.178 

4.199 

26.779 

26.70e. 

BLOCKAGE  FACTOR — OUTLET  TRAVERSE 

■  -0.3142  OUTLET 

RAKE-  ^^3456  ■ 

CALCUL ATEO 

SUDDEN  expansion  AREA 

RATIO— ACTUAL  EXIT 

BLADE  AREA  RATIO  IS 

1.3741 

RADIAL  POSITION 

OUTLET  TRAVERSE 

OUTLET  RAKE 

1 

1.437 

1.299 

2 

0.904 

P.903 

3 

0.927 

0.914 

4 

0.954 

0.954 

_ _5 _ 

0.938 

0.939 

AEDC.TR.69.42 


FRri«<  RUN  C»  31  0094  RFLrCTFO  FX&MPLFS  FOR  TEST  CASES  PAGE  6 

DATA  GROUP  1  l.ON  MIN  EXHIBITING  OUTPUT  FOR  S  STREAMTUBES  WHEN  NN«5 


RESULTS  riF  STREAMTURR  CALCULATIONS 
INLET — OUTLET  T RAVER SP _ 


RADIAL 

total 

AnlABAT IC 

OVERALL 

PASSAGE 

passage 

EXPANSION 

PROFILE 

PROFILE 

DIFFUSION 

PnSITION 

PRESSURE 

EFFICIENCT 

LOSS 

ENTRANCE 

SHOCK 

LOSS 

LOSS 

LOSS 

factor 

RATIO 

MACH  NO. 

LOSS 

PARAMETER 

1 

1.2721 

0.2455 

0.6880 

1.8562 

0.2523 

0.2007 

0.2350 

0.0227 

0.4099 

- 5 - 

1-4419 

0.3546 

0.5621 

1.8318 

0.2415 

0.0180 

0.3026 

0.0350 

0.2967 

3 

1.5429 

0.3942 

0.5550 

1.8060 

0.2310 

0.0082 

0.3158 

0.0381 

0.2891 

4 

1-4239 

0.3419 

0.5843 

1.7776 

0.2194 

0.0370 

0.3279 

0.0375 

0.2993 

5  _ 

1 . 3791 

0.2985 

0.5799 

1.7471 

0.2071 

0.0727 

0.3001 

0.0323 

0.2862 

AV^RAG^ 

1 .40SO8 

0.32693 

0.59334 

1.80373 

0.23024 

0.06732 

0.29628 

0.03312 

0.31627 

INLET  — 

outlet_raxe 

RADIAL 

TOTAL 

ADIABATIC 

OVERALL 

PASSAGE 

PASSAGE 

EXPANSION 

PROFILE 

profile 

DIFFUSION 

POSITIDN 

PRESSURE 

EFFICIENCY 

LOSS 

entrance 

SHOCK 

LOSS 

LOSS 

LOSS 

factor 

RATIO 

MACH  NO. 

LOSS 

parameter 

1 

1.3050 

0.2663 

0.6693 

1,6562 

0.2523 

0.1579 

0.2592 

0.0262 

0.3947 

2  ■ 

1 .4798  " 

0.3761' 

0.5589  ' 

1.8318 

0.2415 

0.0150 

0. 3025 

0.0352 

0.2993 

3 

1.5250 

0.4020 

0.5537 

1.8060 

0.2310 

0.0086 

0.3142 

0.0375 

0.2906 

4 

1.4122 

0.3355 

■  0.'5B57 

1.7776 

0.2194 

0.0407 

0.3256 

0.0370 

0.2997 

5 

1.3297 

0.2972 

0.5801 

1.7471 

0.2071 

0.0724 

0.3006 

0.0324 

0.2862 

AVFRAGF  _ 

1 . 41035 

_  0.  33543 

0,56955 

1.80373 

0.23024 

0.05890 

0-30041 

0-03365 

0.31411 

AEDC.TR. 69-42 


7EST  C4SES  FOR 

STREAMTUBE  PROGRAB 

PAGE  1 

INPUT 

C0NST4H7S 

q 

9.84390 

11.01850 

9.99300 

11.01890  9.90000  11.01850 

A 

0,0 

0^0 

0.0 

0.0  0.0 

0.16605 

0.0 

0.0 

0.14835 

0.0  0.0 

c 

0^90000 

1.00000 

0,0 

0.0 

0.0 

0.0  0.0 

0.05920 

i  “ 

0,0 

0.0 

3.0 

0.01358  0,0 

0.0 

0.0 

6.44387 

-0.50726  0.01734 

FROM  RUN  CA  T1  0094 
DATA  rtnniiP 

SELECTED  EX4HPLES  FOR  TEST  :»SES 
■>  i.nm  MAT  FXHiniTINS  OUTPUT  FOR  10  STRE4HTUBES  WHEN 

NN-11 

CONSTANTS 

FOR  4  SUBSET 

P 

1 1 ,98500 

11.77000 

19.76300 

21.54000  0.0 

0.0 

PBAR 

0-0 

0.0  14. 

9470  14.9430  14.9280 

14.9190 

14.9180 

0,0 

0.0 

0.0 

0.0 

0.9 

0.0 

6.6  2'».T5nO  30.7300 

n.n  0.0  0.0 

31.0030  31. 
0.0  30. 

2500  31.5030  34.0999 
,8930  31.9430  34.5370 

36.2500 

35.9770 

35.9000 
34. 7630 

34.6300 

0.0 

33.5000 

0 

30,1000 

r^o 

0.0 

0.0 

A*  ft 

0.0 

0 

0.0 

TPMP 

776.670 

782.^0 

780.990 

776.010 

771.980 

764.790 

529.790  526.930  526.410 

526.980  528.900  787.670  786.670 

780.670 

774.670 

ANGLF 

Q-»0 

—DEGREES 

0.0  0,0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 _ 

0.0 

0.0 

0.9 

0.0 

90.400  98.190 

0.0  0.0 

86.200 

9.0 

85.000 

85.000 

70.599 

70.500 

55.890 

55.800 

51.000 

51.000 

52.000 

52.000 

52.509 

0.0 

53.100 

0.0 

54.509 

0.0 _ 

0,0 

_o^e_ 

0,0 

0.0 

ft,^ 

HN-16T10. 0000000  00 


AEDC-TR. 69-42 


PAGE 


TEST 

CASES_jaR_  STREAMTUBE 

PROGRAN 

FROM  RUN  C* 

31  DORA 

SELECTED  EXAMPLES  FOR 

TEST  CASES 

DATA 

GR9UP  2 

l.PN  1 

MAX  EXHIBITING  OUTPUT  FDR  10  STREAMTUBES 

WHEN  NN>11 

CALCULATED  BOUNDARIES  AND  CENTERS  OE  STREAMTU9ES 

RATIO  TO  PASSAGE  HEIGHT  AND  RADIUS  IN  INCHES 

J«.ET 

STREAM TUBE 

outlet 

TRAVERSE 

OUTLET  RARE 

STREAMTUBE 

NO. 

OUTER 

CENTER 

OUTER 

CENTER 

OUTER 

CENTER 

1 

INNER 

INNER 

INNER 

1.3000 

3.9?75 

11.0185 

10.9333 

1.0000 

3.7763 

11.0185 

10.7683 

1.0000 

0.7800 

11.0185 

10.7724 

2 

0.8819 

10.8797 

3.6960 

10.6785 

0.7010 

10.6840 

0.8819 

0,8365 

10.8797 

10,8264 

0.6960 

0.6351 

10.6785 

10.6103 

0.7010 

0.64II 

10.6840 

10.6171 

3 

0.7911 

10.7730 

0.5842 

10.5534 

0.5911 

10.5611 

0.7911 

0.7456 

10.7730 

10.7195 

3.5842 

3.5397 

10.5534 

10.5037 

0.5911 

0.5474 

10.5611 

10.5122 

_ 4 _ 

0.6999 

10.6659 

0.4998 

10.4590 

0.5081 

10.4683 

0.6999 

3.6540 

10.6659 

10.6120 

0.4998 

0.4623 

10.4590 

10.4170 

0.5081 

3.4713 

10.4683 

10.4271 

5 

0.6079 

10.5578 

0.4063 

10.3768 

0.4358 

10.3875 

0.6079 

0.5614 

10.5578 

10.5032 

0.4263 

3.3914 

10.3768 

10.3378 

0.4358 

0.4014 

10.3875 

10.3489 

6 

0.5146 

10.4481 

0.3574 

10.2998 

0.3677 

10.3112 

0.5166 

0.6676 

10 .448 1 
10.3927 

0*5574 

3.3240 

10.2998 

10.2624 

0.3677 

0.3345 

10.3112 

10.2741 

7 

0.4198 

10.3368 

3.2910 

10.2254 

0.3016 

10.2373 

0.6190 

0.3719 

10.3368 

10.2805 

3.2910 

3,2581 

10.2254 

10.1887 

P.3016 

0.2668 

10.2373 

10.2006 

e 

0.3236 

10.2238 

3.2252 

10.1519 

0.2359 

10.1639 

0.3236 

3,2750 

10.2238  1 

10.1667  J 

1  0.2252 

0,1921 

10.1519 

10.1149 

0.2359 

0.2028 

10.1639 

10.1269 

9 

0.2261 

10.1091 

0.1587 

10.0775 

0.1693 

10.0894 

0.2261 

•'.1767 

10.1091 

10.0512 

0.1587 

0.1246 

10.0775 

10.0394 

0.1693 

0.1352 

10.0894 

10.0512 

10 

0.1270 

9.9928 

0.3898 

10.0004 

0.1001 

10.0119 

3.1270 

3.3768 

9.9928 

9.9337 

3.0898 

0.0534 

10.0004 

9.9598 

6.1001 

0.0634 

10.0119 

9.9709 

0.0 

9.8435 

0.0 

9,9000 

0.0 

9.9000 

AEDC.TR.69.42 


FROM  RUN  CA  31  00R4  SELFCTED  EXAMPLES  FOR  TEST  CASES  PtGE  3 

DATA  GROUP  2  1.0ft  MAX  EXHIBITING  OUTPUT  FOR  10  STREAMTUBES  WHEN  NN»1I 

SUMMARY  OF  MASS  AVERAGFl  STREAMTU8E  FLOW  VARIABLES  ~ 


INLET 

RADIAL 

POSITION 

_ 1 _ 

RATIO  TO 

PASSAGE 

HFISHT 

O.PZTS 

TOTAL 

PRESSURE 

14.6811 

STATIC 

PRESSURE 

11.7861 

total 

TEMP 

531.5419 

ABSOLUTE 

FLOW  ANGLE 

0,0 

2 

9.B3BS 

14.8924 

11.8050 

529.8412 

0.0 

3 

0.T4S6 

14.9329 

11.8243 

528.5614 

0.0 

4 

0.4549 

14.9495 

11.8436 

527.8224 

0.0 

5 

0.5614 

14.9412 

11.8633 

527.6409 

0.0 

6 

0.4674 

14.9238 

11.8835 

527.7741 

0.0 

T 

r.3Ti9 

14.9174 

11.9038 

528.0092 

9,0 

B 

O.Z750 

14.9145 

11.9244 

528,5386 

0.0 

9  . 

0.1767 

14.9152 

11.9454 

529.3882 

0,0 

10 

0.07SB 

14.0290 

11.9660 

530.5064 

0.0 

AVERAGE 

14.8897 

11.8745 

528.9624 

0.0 

OUTLET 

traverse 

1 

0.77S1 

30.9512 

21.1668 

791.2672 

76.3946 

2 

0.6351 

32.4141 

20.8912 

788.1543 

64,0474 

0.53R7 

33.6540 

20.7213 

785.5254 

57.8011 

4 

0.46?3 

34.8779 

20.5723 

782.5876 

53.9044 

5 

0.3914 

35.7237 

20.4461 

780.5310 

52.0950 

6 

0.3Z4n 

36.1830 

20.3365 

779.5343 

51.1369 

7 

0.Z581 

56. 57?o 

2D.2306 

779,?oqt 

50.7772 

8 

0.1921 

36.3280 

20.1188 

779.7560 

50,9443 

9 

0.1246 

36.0160 

19.9903 

780.9998 

51.6622 

_ 10 

0,0554 

33.4600 

19.8575 

702,8065 

53.0264 

_ 

AVERAGE 

34.5981 

20.4331 

783.0462 

56.2589 

OUTLET 

RAKE 

1 

0.7800 

31.1804 

21,1728 

785.4903 

76.6526 

2 

0.641 1 

32.8269 

20.9025 

783.0721 

65.3194 

3 

0.5474 

34.0746 

20.7352 

781.0166 

58.3364 

4 

0.4713 

35.0766 

20.5901 

779.6713 

54.2292 

5 

0.4014 

35.6826 

20.4633 

778.5305 

52.3195 

6 

0.3345 

35.9406 

20.3533 

777.1745 

51.2513 

7 

0.2688 

35.9334 

20,2478 

775,5259 

50.7067 

8 

0.2028 

35.6723 

20.1374 

773.4523 

50.8575 

9 

0.1352 

35. 1043 

20.0112 

770.7163 

51.5194 

10 

0.0634 

33.518t> 

19.8833 

767,3171 

52.1811 

^  AVERAGE 

34.4810 

2074497 

777.1967 

56.345’ 

AEDC.TR.69.42 


FROM  RUM  CA  31  OOPA  SFtECTEp  EXAMRLE.S  FOR  TEST  CASES _ 

DATA  GROUP  2  l.ON  MAX  EXHIBITIMG  OUTPUT  FOR  10  STREAMTUBES  WHEN  NN«H 

CALCULATED  FLOW  CONDITIONS  ~ 


RADIAL 

. POSITION 

RATIO  TO 

passage 

HEIGHT 

ABSOLUTE 

HACH 

RELATIVE 

MAC*! 

TOTAL 

ENTHALPV 

RELATIVE 

FLOW 

BLADE 

ANGLE 

RELATIVE 

TOTAL 

RELATIVE 

TOTAL 

temperature 

NUMBER 

NUMBER 

ANGLE 

PRESSURE 

INLET 

1 

0.92T5 

0.5619 

1.5629 

127.5659 

68.6532 

61.1529 

47.4402 

7A3.A767 

2 

0.8365 

0.5857 

1.560A 

127.1581 

67.9562 

60.9143 

A7.3AA6 

737.6650 

...  -  3  . 

O.TASi 

0.5871 

I,5A99 

126,85LS_ 

67.7370 

60.6723 

A6.6708 

732.3088 

k 

O.ASAO 

0.586A 

1.5368 

126.67A1 

67.5697 

60.A2A9 

A5.883! 

727.5037 

s 

0,8614 

0.5835 

1,5220 

126.6306 

67.A582 

60,1707 

44.9811 

•'23.2A8A 

6 

n.A6TA 

0.5797 

1.5063 

126.66  2  5 

67.3633 

59.9085 

AA.0359 

719.2850 

T 

0.3719 

0.5769 

1.A90T 

126.7189 

67.230A 

59.6381 

A3.  1263 

715.A062 

8 

0.2T50 

0.5744 

1.ATA7 

126.8458 

67.0773 

59.3591 

A2.2166 

711.8062 

_ _9. 

0.1767 

0,5721 

1,A583 

127,0A95 

66.9016 

_ 59. 071 3 

A1.3079 

708.5116 

10 

0.0768 

0.5620 

1.A3T3 

127.3176 

66.9819 

58.7735 

AO. 1588 

705.A599 

average 

0.5777 

1.5099 

126.9A7A 

67.A929 

60.0086 

AA.3165 

772.A671 

OUTLET 

TRAVERSE 


1 

3.7763 

0.7582 

0.4999 

190.2376 

59.0954 

30.7276 

25.0942 

745.5762 

2 

0.6  351 

0.8 186 

0.57A8 

189,4506 

52.7462 

30.3570 

26.1205 

741.3188 

3 

0.539T 

0.8629 

0.65A5 

188.8115 

45.3744 

30.1053 

27.6136 

742.6115 

4 

0.A673 

0.9030 

0.7072 

188.0974 

41.2108 

29.8997 

28.7023 

740.4232 

5 

0.391A 

O.930A 

0.7330 

187.5976 

38.7638 

29.7111 

29. 209? 

737.1109 

6 

0.3240 

0.9A65 

0.7A58 

187.3555 

37.2144 

29.5309 

29.4015 

734.85"! 

7 

0.2581 

0.9559 

0.7A82 

187.2985 

36.1163 

29.3540 

29.3149 

732.9254 

8 

9.1921 

0.9597 

0. 7409 

187,4093 

35.2999 

29.1764 

28.9531 

731.0246 

9 

0.12A6 

0.9577 

0.7230 

187.7115 

34,7479 

28.9940 

28.2964 

729.2198 

10 

3.0534 

0. 8973 

0.6870 

188.1506 

38.2214 

28.8011 

27.2126 

738.1598 

AVERAGE 

0.8993 

0.681A 

188.2090 

4?, 6790 

29.6657 

27.9918 

7^7.^734 

OUTLET 

RAKE 

1 

0.7803 

0.7656 

P.A975 

IBB. 8029 

69.1890 

30.7372 

25.0620 

738.2983 

2 

0.6411 

0  .  8305 

0.5685 

188.2152 

52.4125 

30.3730 

26.0133 

733.021? 

3 

0.5A74 

0.87*0 

0*6466 

187.7156 

44.9840 

30.1255 

27.4962 

734.8428 

4 

3.4713 

0.907A 

0.70AA 

187.3887 

41.1491 

29.9237 

28.6569 

736.1385 

5 

O.A0I4 

0.9286 

0.7310 

187.1116 

39.0609 

29.7376 

29.1785 

735.2389 

6 

0.33A5 

0.9399 

0.7AA9 

185.7821 

39.8436 

29,5588 

29.4013 

734.0343 

7 

0-2688 

0.9AA3 

0.7A88 

186.3817 

37.1264 

29.3827 

29.3551 

732.1858 

9 

3.2028 

0.9A27 

0.7A32 

185.8781 

36.8010 

29.205? 

29.0418 

729.5112 

9 

0.1352 

0.9339 

0.7270 

185.2137 

36.9336 

29.0226 

76.471P 

725.8543 

10 

0.063A 

0.8921 

0.7077 

184.3885 

39.3780 

28.8281 

27.7554 

728.4797 

AVERAGE 

0.9959 

0.6822 

186.7878 

43.4878 

29.6895 

28.0391 

732.7605 

••‘CE  _  A 
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comnfnts 


PAGF 


HASS  FLOW 


PERCFNT  HASS  LOSS  IN  OPTFRHIHtNG  RADIAL  POSITIONS 

PFRCFNT  ENERGY  UNBALAN:E 


INLET 

?O.Tl ?D  no 


?.lYno-OB 


OUTLET  TRAV. 


?0,6A«>0  00 


STRFAHTUBF 


OUTLET  TRAV 


T. 0850-07 


OUTLET  RAKF 


OUTLET  RAKF 
20. 3870  no 


3.0080-07 


6 

7 

B 

<> 

ID 


DELTAIHI 

67.542 

52.202 

61.050 

61.423 

60.057 


60.503 

60.580 

60.554 

50.562 

60.833 


DELTAIEI/OELTAIHI 
3.730120  00 
5.074860  00 
7.648990  "0 
8.402440  00 
_ 8.371120  nn 


8.550550  00 
8.757850  00 
8. 851720  00 
8.857540  00 
1.303230  01 


OELTAIHI 
61.237 
51.057 
60.854 
'  50.715 
50.481 


60.120 

59.653 

59.032 

58.154 

57.071 


OFLTAIFI/OFLTAIHI 
_8. 848910-01 
1.796870  no 
_  4.574170  on 

6.544380  on 
_ 7.525680  00 


8.124580  00 
8.413B10  00 
8.770210  on 
7.53554D  no 
9.501060  00 


average 


61. 262 


8.128560  on 


59.840 


6.346510  no 


STATIC  PRESSURE  AT  TRAILING  EOGE 


U1 

03 


EROH  NORHAL  SHOCK  EQUATIONS  FROM  SUOOEN  EXPANSION  EQUATIONS 


BFLATIVF  total  pressure  at  trailing  EOGE 
FROM  SUOOEN  EXPANSION  EQUATIONS 


STREAHTUBE 

1 

2 

3 

4 

5 


6 

7 

8 
9 

10 


PASSAGE  entrance 
44.  807 
44.311 

43.828 

43.332 

_ 47.735 _ 

47.142 

41.536 

40.905 

40.250 

39.606 


OUTLET  TRAV. 
17.965 
15.493 
16.147 
16.317 
16.645 
15.921 
17.159 

17.363 
17.509 

16.363 


OUTLET  RAKE 
18.009 
16.630 
15.330 
16.395 
_ 16.588 


16.712 

15.774 

16.766 

16.559 

15.798 


outlet  TRAV. 
_  32.884 

29.535 
29.546 
29.8*4 
29.963 


29.947 

29.734 

29.289 

28.582 

27.964 


OUTLET  rake 
32.785 
29.275 
29.281 
29.782 
29.998 


30.041 

29.914 

29.574 

28.995 

28.764 


BLOCKAGE  FACTOR— OUTLET  TRAVERSE-  0.1843  OUTLET  RAKE> - 0.1891 - 

CTlCULATEO  SUOOEN  EXPANSION  AREA  RATIO— ACTUAL  EXIT  BLAOE  AREA  RATI 0  IS 


1.3741 


RAOIAL  POSITION 
_ 1 


OUTLET  TRAVERSE 
_ 4.229 


2 

3 

4 

5 

6 
7 

_  o 

io 


1.963 

1.490 

1.305 

1.223 

1.181 

1.158 

1.145 

1.142 

1.257 


outlet  rake 

4.152 _ 

1.945 


1.478 

1.304 

1.232 

1.197 

1.181 

1.180 

1.195 

1.278 
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FROM  RUN  CA  31  0094  SFLFCTEO  FXAMPLFS  FTP 
DATA  GROUP  2  1,0N  MAX  CXHIRITING  OUTPUT  i 

RFSULTS  OF  STREAMTUBE  CALCULATIONS 

INLET — OUTLET  TRAVERSF 

TFST  CASES 

FOR  10  stpsahtubes  when 

NN»11 

RADIAL 

total 

adiabatic 

OVERALL 

PASSAGE 

PASSAGE 

FXPANSION 

oRDFIlF 

PROFILF 

DIFFUSION 

POSITION 

PRESSURE 

EFE  ICIENCY 

LOSS 

ENTRANCE 

SHOCK 

LOSS 

LOSS 

LOSS 

factor 

RATIO 

MACH  NO. 

LOSS 

PARAMETER 

I 

2.1082 

0.4814 

0.5873 

1.0442 

0.2479 

0.2185 

0.1209 

0.0073 

0.7152 

2 

2.1766 

0.5057 

0.5459 

1.8330 

0.2421 

0.0961 

0.2077 

0.0211 

0.6605 

3 

2.2537 

0.5327 

0.4956 

1.8219 

0.2375 

0.0555 

0.2026 

0.0237 

0.5981 

4 

2.3331 

0.5625 

0.4581 

1.8100 

0.2325 

0.0338 

0.1917 

0.0239 

0.5570 

5 

2.3910 

0.5851 

0.4364 

1  .7972 

0.2274 

0.0228 

0.1863 

0.0240 

0.5353 

6 

2.A245 

0.5986 

0.4236 

1.7838 

0.2220 

0.0170 

0.1847 

0.0241 

0.5219 

7 

2.4383 

0.6044 

0.4199 

1.7701 

0.2164 

0.0134 

0.1901 

0.0251 

0.5155 

R 

2.4358 

0.6043 

0.4251 

1  .tSSB 

0.2105 

0.0111 

0.2035 

0.0270 

0.5162 

9 

2.4147 

0.5976 

0.4402 

1.7409 

0.2044 

0.0097 

0.2261 

0.0300 

0.5246 

10 

2.2564 

0.5456 

0.4657 

1.7261 

0. 1989 

0.0266 

0.2402 

0.0303 

0.5373 

AVfRAGE 

2.32321 

0.56180 

0.46980 

1.78829 

0.22396 

0.05045 

0.19539 

0.02364 

0.56817 

INLET--OUTLET  RAKE 

radial 

TOTAL 

ADIABATIC 

OVERALL 

PASSAGE 

PASSAGE 

EXPANSION 

PROFILE 

PROFILE 

01 FFUSION 

POSITION 

PRESSURE 

efficiency 

LOSS 

ENTRANCE 

SHOCK 

LOSS 

LOSS 

LOSS 

factor 

RATIO 

MACH  NO. 

LOSS 

PARAMETER 

1 

2.1230 

0.49e0 

0.5892 

1.8442 

0.2479 

0.2166 

0.1247 

0.0075 

0.7194 

2 

2.2043 

0.5254 

0.5505 

1.8330 

0.2421 

0.0904 

0.2180 

0.0223 

0.6685 

3 

2.2ftl9 

0.5517 

0.5009 

1.821° 

0.2375 

0.0512 

0.2122 

0.0250 

0.6060 

4 

2.3463 

0.5734 

0.4618 

1.8100 

0.2325 

0.0331 

0.1962 

0.0245 

0.5609 

5 

2.3882 

0.5890 

0.4400 

1.7972 

0.2274 

0.0242 

0.1885 

0.0242 

0.5373 

6 

2.4083 

0.5902 

0.4?64 

1.783B 

0-2220 

0.0199 

0.1845 

0.0239 

0.5221 

7 

2.4088 

0.6042 

0.4215 

1.7701 

0.2164 

0.0179 

0.1873 

0.0244 

0.5142 

B 

2.3918 

0.6058 

0.4251 

1.7558 

0.2105 

0.0176 

0.1970 

0.0256 

0.5132 

9 

2.3536 

0.6030 

0.4386 

1  .7409 

0.2044 

0.0192 

0.2150 

0.0277 

0.5201 

_ 10 

2.2468 

0.5784 

0.44^3 

1.7261 

0.1989 

0.0358 

0.2146 

0.0266 

0.5233 

AVERAGE 

2.31538 

0.57283 

0.47033 

1.78829 

0.22396 

0.05257 

0.19380 

0.02317 

0.56850 
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TEST  C»SES  EfIR  STREARTUBF  PROGRAM 


PAGE 


1 


INPUT  constants 


R  ■ 


'll  .01  950 

9.90000 

11.01850  9.90030 

11.01950 

A 

C.O 

_ 6 

0.0 

0.0 

0.0  0.0 

0,16605 

0.0 

0.0 

0.16935 

0.0  0.0 

c 

O.OOPOO 

1. 00030 

0 

0.0 

0.0 

0.0 

0.0  0.0 

3.05520 

E 

P.P 

E 

0.3 

0,0 

0.01358  0.0 

O.P 

0.3 

6,44^fl7 

-0.50725  0.01736 

- 

FROM  RUN 

CA  31  0096 

SELECTED 

EXAMPLES  EOR  TEST  CASES 

DATA  GROUP 

3  l.ON  MAX  EXHIBITING  OUTPUT 

FOR  5  STREAMTUBES  WHEN 

NN-13 

CONSTANTS 

FOR  A  SUBSET 

P 

ll.P90''0 

11.75330 

19.67331 

19.00565  19.78331 

21,56030 

PflAR 

0.0  0.0 

0.3 

3.0  0. 

0  16.8550  16.9633 

16,9267 

16.9150 

16.9133 

0.0 

0.0 

f>.0 

O.ft 

0,0 

0.0  ?5.3000  27.5300 

28.2500  28. 

5003  28.5000  29.2000 

36,7030 

37.8500 

36.0000 

35.5000 

114. POO? 

36.1033 

32.55np 

3.3 

0.0  0.0 

0.0 

0.0  0. 

0  3O.92P0  31.9833 

36.5766 

36.0100 

36.8066 

O.P 

0.0 

0,0 

0.3 

3.0 

TFMP 

529.736  526.690  525.271 

527.051  529 

.265  787.670  785.670 

780.670 

>76.670 

776.670 

782.728 

781.055 

775.976 

771.852 

755.216 

ANGLE  — OEGREES 

0.0  0.0 

0.0 

O.P  0 

•  0  0,0  o»n 

P.O 

0.0 

0.0 

0.0 

0.0 

0.0 

P.O 

3.0 

0.0  108.600  106.630 

n?.47n  99 

.730  99.600  69.600 

67. 730 

67.900 

51.300 

51.900 

52.503 

53.303 

3.0 

0,0  0.0 

0.0 

0.0  0 

.0  85.0P0  70,503 

55.800 

51.000 

52.000 

0.0 

0.0 

O.P 

O.P 

0.0 

RN-16727.: 

3300000  00 

AEDC-TR.69.42 


c;i 


^127?  _ 9,9930  0,09»7 _ 9^99^  _ 0,10^  10,0120 

0.0  9,8435  0.0  9,9000  0.0  9.9000 


AEDC-TR. 69-42 


P»G= 


FROM  RUN  C*  3J_  009*  SELFCTED  EXRMPtES  FOR  TEST  OSES _ 

DAT*  GROUP  D  l.ON  MAX  EXHIBITING  OUTPUT  FOR  5  STREAHTUBES  WHEN  NN>13 

iuMMARr'WMASS  Tvf“rAGEO~STr'eAMTUBE  “FLON  VAIMABLES 


T 


RADIAL 

RATIO  TO 

TOTAL 

STATIC 

TOTAL 

ABSOLUTE 

POSITION 

PASSAGE 

PRiSSURE 

..PRESSURE 

TEMP 

FLOW  ANGLE 

HEIGHT 

'  INLET 

1 

0.5820 

1*.T923 

11.790* 

530.6331 

0.0 

2 

0.7002 

l*.9*l* 

11.8313 

528.  0566 

0.0 

3 

J5.515S 

I*.93l5 

LL.8732 

527.58*7 

0.0 

4 

0.32*0 

1*.9158 

11.9166 

528.3295 

0.0 

5 

J  JiTi 

1*.86B9 

11.9609 

530.27J5 

0.0 

AVERAGE 

1*.8900 

11.87*5 

528.9753 

0.0 

OUTLET 

TRAVERSE 

_1 _ 

0.5997 

35.3261 

18.7550 

787,6370 

**.1672 

2 

0.*678 

37.07*5 

18.7288 

785.2097 

*7.1127 

3 

0.3531 

37.7026 

18.6667 

781. *266 

*7.2920 

4 

0.2236 

37.1536 

18.5922 

780.3389 

*9. *239 

5 

0.08*7 

35.3806 

18.5180 

782.7836 

51.888* 

_  AVERAGE 

36.527* 

18.6521 

783.5792 

*7.976R_ 

OUTtFT 

RAKE 


_JL  . 

0.7010 

32.0*21 

21.0*35 

78*. *0*7 

70.9886 

2 

0.5082 

3*. 6129 

20.6728 

780.2107 

56.2870 

3 

0.3677 

35. 8**7 

20. *213 

777.6833 

51.7866 

4 

0.2360  • 

35.8381 

20.2081 

77*. 5021 

50.8219 

_ 5 _ 

0.100! 

3*. 2526 

19.9655 

769.  *693 

51.8500 

AVERAGE 


3A.51R1 


20.JA622 


777.25*0 


56.3*68 


AEDC-TR. 69-42 


FROM  RUN  CA  00R4  SFIFCTFO  FXAMPLES  FOR  TEST  CASES 
HATA  group  3  l.ON  HAX  EXHIBITING  OUTPUT  FOR  5  STREANTUBES  WHEN  NN*13 

CALCULATED  FLOW  CONDITIONS 


RADIAL 

RATIO  TO 

ABSOLUTE 

RELATIVE 

TOTAL 

RELATIVE 

BLADE 

RELATIVE 

RELATIVE 

POSITION 

PASSAGE 

NACH 

NACH 

ENTHALPY 

FLOW 

ANGLE 

TOTAL 

TOTAL 

HFIGHT 

NUNBER 

NUNBER 

ANGLE 

PRESSURE 

TFMPFPATIJPF 

INLET 

— 

— 

1 

3.W23 

0.5784 

1.5637 

127.  3480 

68.2897 

61.0342 

47.5144 

740.9455 

? 

0.7D02 

0.5870 

1.5449 

125.7302 

67.66R1 

60.5503 

46.3845 

75^.7055 

3 

3.5153 

0.3815 

1.5157 

126.6171 

67.45T5 

60*041 n 

44*6044 

721 .5656 

4 

3. 3240 

0.5754 

1.4838 

126.7957 

67.1852 

59.5007 

47.7491 

714.0634 

5 

0.1272 

0.5662 

1.4483 

127.2615 

66.9865 

58.9245 

40,7759 

707.6908 

AVERAGE 

0.5777 

1.5113 

126.9505 

67.5135 

60.0103 

44.4057 

722.8041 

OUTLET 

TRAVERSE 

1 

0.5997 

0.0964 

0.8924 

189.3248 

36.7854 

30.2637 

31.4463 

761  .9600 

2 

3.4675 

1.0385 

0.8438 

188.7347 

33.1181 

79.9144 

29.  8293 

73fi.''954 

3 

0.3501 

1.0552 

0.8386 

187.9368 

31.4138 

29.6007 

29.5720 

7?9.6‘»47 

4 

0.2236 

1.0463 

0.7897 

187.5510 

30.4732 

29.2613 

28.0398 

720.2505 

5 

0.0847 

1 . 0086 

0.7262 

188.1450 

31.0075 

28.8859 

26. 2919 

719.3741 

AVERAGE 

1 . 0290 

0.8182 

188.3385 

32.5588 

29.5852 

?9.054T 

733.8744 

nUTLET 

RAKE 

1 

0.7013 

0.7998 

0.5213 

188.5391 

60.0146 

30.5304 

25.3147 

755.6599 

2 

0.5082 

0.8914 

0.6762 

157*5198 

47.9679 

33.0216 

7A.0645 

755.0T16 

3 

0.3677 

0.9345 

0.7387 

186.9057 

38.5073 

29.6478 

29.3284 

734.5577 

A 

0.2363 

0.9437 

0.7467 

186.1330 

37.0272 

29.2945 

20.2416 

750.964! 

5 

0.1001 

0.9138 

0.7177 

184.9109 

38.1447 

78.9277 

?A.l?5? 

T27.525T 

AVERAGE 

0.8966 

0.6801 

186.8017 

43.3323 

29.6844 _ 

28.0149 

732.3558 

PaG=  * 


AEDC.TR.69.42 


C0MWFNT5 


P»GE 


5 


_  _  _ _ _ IW-PT 

mss  FLOW _ ?O.TnD  00 

PERCENT  NASS  LOSS  IN  OETERNINING  RADIAL  P0SITI3NS  1.047D-08 

PERCENT  ENERGY  UNBALANCE 


OUTLET  TRAV.  OUTLET  RAK= 

?I.fe8»n  00 _ PO.AOTO  no 

A.985D-P7  UB^Pn-OT 


STREANTUBE  OUTLET  TRAV 


OUTLET  RAKE 


CJI 

CD 


2 

3. 

4 

5 


OELTAIMI  DELTAIEI/DELTAIHI 

61.R77  _  Il327460  01  _ 

62.004  7.117J30  00 


61.320 _ 5.968970  00 

60.755  4.012160  00 

60.884 _ 5.072550  00 


OELTAIMI 

61.191 

60.790 

60.289 

59.337 

57.649 


OELTAIEI/DELTAIHI 
5.679540-01 
5.357070  00 

7.691270  no 

8.213290  90 
_ 8.454050  00 


_  AVERAGE  61.388 _ 7^088520  00  _  59.851  6.056770  00 

_ static  PRESSUREJAJLTRAILING.EDGE  RELATIVE  TOTAL  PRESSURE  AT  TRAILING  EDGE 

FROM  NORMAL  SHOCK  EQUATIONS  FROM  SUDDEN  EXPANSION  EQUATIONS _ FROM  SUDDEN  EXPANSION  EQUATIONS 


STREANTUBE 

PASSAGE  ENTRANCE 

OUTLET  TRAV. 

OUTiET  RAKE 

OUTLET  TRAV. 

OUTLET  RAKE 

1 

44.634 

14.986 

17.175 

31.843 

30.271 

2 

43.661 

16.979 

16.345 

29.920 

79.494 

3 

42.540 

17.664 

16.627 

29.600 

30.058 

4 

41.314 

17.988 

16.750 

28.052 

29.798 

5 

40.008 

17.623 

16.222 

26.326 

79.894 

BLOCKAGE  FACTOR— OUTLET  TRAVERSE-  -6.2874  OUTLET  RAKE-  0.1492 


CALCULATED  SUOOEN  EXPANSION  AREA  RATIO— ACTUAL  EXIT  BLAOE  AREA  RATIO  IS  1.3741 


RADIAL  POSITION  OUTLET  TRAVERSE  OUTLET  RAKE 
_ 1  _ L.095  _  2.673 

2  1.058  *  1.382 

3  _  1.034  _  1.215  _ 

4  1.026  1.182 

_ 5 _ 1.032 _ 1.23* _ 
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FROM  RUN  C*  31  OOR4  SFLFCTFO  FXANPLFS  FOR  TFST  CASES 
DATA  GROUP  3  l.ON  MAX  EXHIBITING  OUTPUT  FOR  5  STREAMTUBES  WHEN  NN=13 


RESULTS  OF  STREAMTU8E  CALCULATIONS 


INLFT— OUTLET  TRAVERSE 


'  P^OIAL 

TOTAL 

adiabatic 

overall 

PASSAGE 

PASSAGE 

EXPANSION 

PROFILE 

PROFILE 

DIFFUSION 

POSITION 

PRESSURE 

EFFIClENCr 

LOSS 

ENTRANCE 

SHOCK 

LOSS 

LOSS 

LOSS 

FACTOR 

RATIO 

MACH  NO. 

LOSS 

PARAMETER 

1 

2.3Bfll 

0.5779 

0.3772 

I.fl404 

0,2458 

0.0113 

0.1202 

0.0162 

0.4331 

2 

2,4813 

0.6036 

0.4213 

i.eisi 

0.2359 

0.0026 

0.1828 

0.0255 

0.4717 

3 

2.5250 

0.6233 

0,4213 

1.7926 

0,2256 

0.0009 

0.1949 

0.0274 

0.4680 

4 

■?.490R 

0.6194 

0.4589 

1.7647 

0.2143 

0.0004 

0.2442 

0.0343 

0.4950 

5 

_ 2-3T9_5 

0.5 85 2 _ 

0.5031 

1.7348 

0.2023 

0.0012 

0.2996 

0.0413 

0.5278 

AVER  AO? 

2.4529P 

0.60187 

0.43638 

1.79010 

0.22477 

0^00327 

0.20834 

0.02893 

0.47913 

INLET— 1 

nuTLFT  RAKE. 

' 

RAOIAL 

TOTAL 

ADIABATIC 

OVERALL 

PASSAGE 

PASSAGE 

EXPANSION 

PROFILE 

profile 

DIFFUSION 

POSITION 

PRPSSORF 

EFFiClENCY 

LOSS 

ENTRANCE 

SHOCK 

LOSS 

LOSS 

LOSS 

FACTOR 

RATIO 

MACH  NO. 

LOSS 

PARAMETER 

1 

2,1661 

0.5121 

0.5748 

1.8404 

0.2458 

0.1387 

0.1903 

0.0150 

0.7040 

2 

2.3165 

0.5633 

6.482‘9 

■  1.8181 

0.2359 

0.0414 

0.2056 

0.0250 

0.5844 

2.4006 

0.5948 

0.4335 

1.7926 

0.2256 

0.0223 

0.1856 

0.0239 

0.5297 

A 

2,4027 

0.6060 

0.4233 

1.7547 

0.2143 

0.0180 

0.1910 

0.0248 

0.5135 

5 

2,3036 

0,5922 

0.4438 

1.7348 

0.2023 

0.0267 

0.2148 

0.0272 

0.5217 

AVERAGE 

2.317P3 

0.57368 

0.47166 

1.79010  _ 

_ 0.22_477 

0.04942 

0.19746 

0.02339 

0.57068 
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